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Re „Fir ah Fal, daß die er Ber M eßpunkte im n-dimensionalen. Pride in. Wäse % : 

| n einer unbekannten Anzahl unabhängiger Variabler mit Gaußischer Verteilung abhängen, jedoch 
treuungsiensor aus Versuchen bekannt ist, wird der von den Meßpunklen. erfüllte Teilraum da di 
cheinlichkei sdichte in demselben sowie die Genauigkeitsh hyperfläche bestimmt. i 


3 In ihe case that the coordinates of the measuring. points (in a n-dimensional space)-are depending Beni 

 — in an unknown way — on an unknown number of independent variables distributed according to the 
G@auss law, that ihe lensor of variation, however, is found out by.ewperiment, the part of Ihe, space filled 
x with Berg | ne ER er density u ‚probability in that part as well as the SurensVr[aon, of er r 
‚are eva unle i 5. 
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AR HENSBECTHEIM o6pasom BABHCHT OT  HEHSBECTHOTO yncJa HE3ABUCHMEIX TIEPeMEHHHIX, PAac- 
.  DpeneleHHBIX 0 .3akoHy T'aycca, Ho Korna TeHs0p paccegauna W3BeCTeH M3 OMEITOB, 
za en YACTb HPOCTPAHCTBA, BAIONHCHHOTO TOYKAMH H3MEPEeHHA, TYCTOTA ee 3 
 B 3TOM ee a en TuNEePNOBePXHOCTb TONHOCTH AsMepenun. a 


$1. Problemstellung in DT 


ai kommt in den Anwendungen oft vor, daß das Ergebnis einer Messung aus n Werten 

En 2 5 %, also einem Vektor = (x, ..., %,) im n-dimensionalen Raume, besteht derart, 
daß die dy Funktionen von r aleatorischen Variablen «,, ..., u, sind, die unabhängig von- 

einander je einer Gaußischen Verteilung mit dem Erwartungswerte Null?) genügen.‘ Die genaue. ie: 
Gestalt dieser Abhängigkeit ist dabei unbekannt; dagegen kennt man aus Versuchen den sog. R 
Streuungstensor T von x, dessen Komponenten w;, durch die Erwartungswerte E(x;x%;) der n 
Produkte &;x; gebildet werden. Oft genügt es, an Stelle der genauen Abhängigkeit der x, von 

- den «, eine lineare Näherung zu verwenden, so daß vorausgesetzt wird, daß zwischen den Vek- 
toren r und u = (w, ..., w,) ein affiner Zusammenhang besteht. Damit ergibt sich die folgende 

Se Aufgabe: Es sei bekannt, daß 


EM N een (1) 


ist, wo X eine unbekannte Matrix mit n Zeilen und r Spalten bei unbekanntem r ist und die 
Komponenten von u unabhängig voneinander je einer Gaußischen Verteilung gehorchen. Bekannt 
sei dagegen der Streuungstensor I. Die Wahrscheinlichkeitsdichte op(x) und die Gleichung der 
Genauigkeitshyperfläche im g-Raum sind zu finden. FA 
Diese Aufgabe wurde kürzlich von K. Stange?) in einer ausführlichen Arbeit für n=3 Se 
behandelt. Es sind dabei außer der unwesentlichen Erweiterung auf beliebiges n noch die folgen- z 
den Fragen offen ‚geblieben: . | % 
a) Es wurde vorausgesetzt, daß T nicht singulär ist, d. h. daß seine Determinante |T| nicht : 
verschwindet. Geometrisch bedeutet dies, daß die g den gesamten r- -Raum erfüllen, was auch 
bei r>n nicht der Fall zu sein braucht. Der Fall |T|=0 muß also noch behandelt werden. 
b) In den Fällen r<n wurde angenommen, daß A und damit der durch die x aufgespannte 
Teilraum R bekannt ist, so daß man in demselben neue Koordinaten einführen und damit n 
"reduzieren kann. Es bleibt also die Frage der Bestimmung von R und g aus T bei r<n offen. 
c) Es ist bemerkenswert, daß im Unterschied zu der;von Stange für 9 angegebenen Formel 
seine Gleichung für die Genauigkeitshyperfläche auch bei T]=0 und r<n richtig bleibt, ob- 
gleich in seinem Beweis |T| +0 und r> n verwendet wird. Der Beweis dieser allgemeinen Gültig- 
keit ist also noch nachzuholen. 
Es erscheint daher zweckmäßig, die Fragestellung nochmals äufzugreifen. 


1) Diese Normierung ist natürlich keine Einschränkung der Allgemeinheit. 
2) K.Stange, Über die Verteilungsdichte der Meß- oder Beobachtungsfehler eines dreidimensionalen 
Punktraumes. Z. angew. Math. Mech. 23 (1948), S. 235—243. ; 
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Nun ist, wenn der Querstrich über einer Matrix die Vertauschung von Zeilen und Spalten bedeutet 


en. 


Es ist Habe Einschränkung der Allgemeinheit, wenn wir im folgendeı 
daß die w, die Streuung E(w)=1 besitzen. Dies läßt sich ja ‘ohne Änderı sch 
lichkeitsverteilung der z.durch Multiplikation der «, mit einem Faktor und gleichzeitige Div 
der entsprechenden Spalte von A durch denselben erreichen. (1) bedeutet nun geometri | 
- der r-dimensionale u-Raum affin auf einen VE IRRSMURGER UntarTene R des gRaumes 12 
m<n und m <r abgebildet wird. ee | 
Betrachten wir zunächst den einfachsten Fall m=n—=r und aka ||+0. Dasich de 
Volumeneinheit bei der affinen Abbildung mit |X| multipliziert, gilt hier für die WpraREIgiehr er 
keitsdichten p(x) und y(u): 


e Bag 


WM) = y (u) = 8 N‘ R 


ww Ar Aı,=- 2: PYr-r-(AA)r. 

Andererseits ist 

E(%)= E (20% Ak Ug Ye) — PALFPLTPR 
0,0 o 


also g 
T=-AA und | = MENT ETTAOE AE (2). 
Damit wird schließlich « 


3 5 
Kerr in. © 3). 
In diesem einfachsten Falle ist die Aufgabe a der Wahrscheinlichkeitsdichte also bereits 
gelöst. 
Bemerkung 1: Ist X speziell eine orthogonale Transformation, so wird T= AA= € 
(= Einheitsmatrix). (3) zeigt dann, daß bei orthogonalen Transformationen wieder unabhängige 
normierte Gaußische Variable entstehen. 

Gehen wir nun wieder zum allgemeinen Falle über. Hier wird der durch Au=0 definierte 
(r—m)-dimensionale Unterraum U des u-Raumes bei der Abbildung annulliert. Wir wählen 
eine orthogonale Transformation R des u-Raumes, po=XRu, bei der U geradein den durch die 
letzten »—m Einheitsvektoren aufgespannten Teilraum übergeführt wird. Es ist dann: 


= Au= ARD = Br, 


wo p nach Bemerkung 1 die gleichen Voraussetzungen wie u erfüllt und in ® die letzten r— m i 
Spalten verschwinden. Da somit YO41s-.., % gar Er t beeinflussen, können wir sie weg- ; 
lassen und statt des v-Raumes gleich «den d turch Vs + % definierten m-dimensionalen Raum ei 
bei m<n verwenden. Wir haben somit den 


„Satz 1: Unter den in der Aufgabe gestellten Voraussetzungen ist = 3b, wo die Kom- 
ponenten ©, ..., m unabgängig voneinander Gaußischen Verteilungen mit der Streuung 1 
gehorchen und ® den Rang m besitzt. Es ist m<n. 


$ 3. Die Wahrscheinliehkeitsdiehte 

Für mn ist |®8|++0 und damit nach (2) auch |T|= |8]?+0. Für m<n dagegen ist 
IT|= |8%]= 0. Im nichtsingulären Fall |T|=+0 ist also sicher m=n, so daß (3) anwendbar 
ist. Wir haben somit 

Satz 2: Bei |T|=#0 ist 

a . . 
Ye Ne in 
Y2r'yiz 

Damit haben wir das Ergebnis von Stange auf etwas einfachere Art wiedergefunden. 


Wir wenden uns nun zu dem offen gebliebenen Falle |T]=0, resp. m <n (Fragen (a) und (b) 
von $1). Hier wird vermöge t=®Bv der m-dimensionale v-Raum auf den m-dimensionalen 
Unterraum R des g-Raumes abgebildet. Wir wählen nun eine‘orthogonale Transformation Ns 


BR 


IR Y N Ur uns 


a ‚Richter, Zur Gau ßischen Verteilung im n-dimensionalen Raume 163 


die Rin den durch die ersten m-Einheitsvektoren aufgespannten Raum R, dreht. Es ist dann bei 


Y= MB < Co 


(6) an 


wo &, eine m-reihige quadratische Matrix mit |E,|#0 ist. 
Ink, wird damit y repräsentiert durch den m-dimensionalen Vektör = (Y1 +. Ym) = Ep 
mit der | durch Satz 2 gegebenen m-dimensionalen Verteilungsdichte | 


mit 


(9) h ZU % 0 
I ee 
: 2m ET 


Da %, orthogonal ist, ist diese Wahrscheinlichkeitsdiehte gleich der entsprechenden Wahr- 
scheinlichkeitsdichte im. Unterraum R des r-Raumes: 


bei 0,8, und. 002761 


| DIENEN AU SEEN a Ne (7). 
Nun ist für yEk.: * 
N | 30 
Hy = y° (% 0) 2 nur: (9 Or, 
oder 
| ee 0 
ty r:H(d 0) Mar. RE N pe N ee (8). 
Andererseits ist nach (5): 
T- BE —- NCCH,-R, (2 0) DE N (9). 


Wird umgekehrt T durch eine beliebige orthogonale Transformation €, in die Gestalt gebracht: 


PBE ca Ki 6) En Tol+0, 


so ist wegen der positiven Eigenwerte von %o: 


= or ) SIERT Ne und. le ee (10), 


wo R.„ und R,-m Tesp. m- und (na—m)-dimensionale orthogonale Transformationen sind. Es ist 
dann nach (8) a 


-1 

RE su ©, Fo: s u REN RL EL AU (11). 
NR 0 
0 Rum 
m-Einheitsvektoren. Da aber der rechts stehende Faktor dieses Ausdruckes für R, den Raum Ir 
in sich dreht, genügt es statt dessen, ©, auf R, anzuwenden. Diese Tatsache liefert zusammen 
mit (6), (7), (10) und (11) den 

Satz 3: Bei |{|=0 bringe man X durch eine orthogonale Transformation ©, auf die Gestalt 


2-8, %0)e, 


wo T, vom Grade m mit X|=+ 0.ist; dann liegen alle y in dem durch die ersten m Spalten von ©, 
aufgespannten Teilraum R und haben dort die m-dimensionale Wahrscheinlichkeitsdichte 


Der Unterraum R der r entsteht durch Anwendung von R,—= € ( auf die ersten 


. 1 
1 SET — (T710 
m m .* 2 bei % = C, ( 8 )e. Q 
a y2 ne" Yz : ih, 


Bemerkung 2: Setzt man in die Formel von Satz 3 einen beliebigen Vektor ein, so liefert 
op(x) die Wahrscheinlichkeitsdichte in dem‘Punkte der Projektion von & auf R. 
IT® 
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„san in belhe Vkar m n y-R 
ln de PROBINDE. DON TERE: ' da 
' genommen lea, die Wahrscheir it 
“ "nktensyetamn, bei welchem n der erste. svekti 
0 Ist n senkrecht zum Unterraum AR, so ist natürlich gu = 0. 
Projektion in R, so setzen wir analog zu (4): Be ER h H 
‚die Hyperebene 4 - BE ee un. Unterraum R, die I 
' gemeinsam hat, wo m;+0 die Projektion von m in R, bedeutet. Der eb 
e entspricht N al Hyperebene KCeth: SR Cm, die Dir d 
ER | ET VE ER LE 
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vom natur besitzt. Da im v-Kaume orthogonale RR nichts 
machen (vgl. Bemerkung 1), können wir annehmen, daß der Normalenvektor : Sm Ka 
Grundvektor ist. Wegen der Affinität aller genannten PETARTER, ist also die gesuchte W 


‚scheinlichkeitsverteilung al F 

BER Ba am P a 
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2 EN \ . ; ern, ©,C,m J my: GE, m ' A j 
Ar 58 Diese Formel haben wir nun noch durch einen Ausdruck in T/zu erseizen. Nun ist Gemäß (6) 
2 a und (9) für beliebiges n: | ig 
e my,’ E,&,m, = m; ae ee = m: (6 5 0)" Rn (6° 0) Kunz N E x | 
us oder | a 

Da 2: % > 
m’ &,&,m, =tn: Rn 00 Run FAR" In Nee BEN us 


I aindere ist n-Tn—=0 für m — 0, also für die n senkrecht zu R, für welche RT, jet, 
Damit haben wir den 


Satz 4: Ist nein beliebiger Vektor der Länge 1, so ist die Wehrseheipiichpeiksäinite von rn: 


1 
m ea "1.6 . & 
Ir y2 T Yn h T n 0 
Yn ist eine Gaußische Verteilung mit dem Genauigkeitsmaß | 


# 


1 | 
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wu 


Damit ergibt sich sofort der’ 
Satz 5: Trägt man in jeder Richtung n das zu Qu gehörige Genauigkeitsmaß Ah. ab, so liegen - 
die Vektoren 3= hyn:n auf der Hyperfläche 
1 
8. Te—= D% | 
Damit ist auch die Gleichung der Genauigkeitshyperfläche als allgemein gültig bewiesen. 
Eingegangen: 26.9. 1948. 
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BORe he valuation of the grors given by T ollmien for the degree of approximation ey 3 in the case of inter- 
RE > _polation and extrapolation, is generalized to any degree of approximation and made more distinct, too. Re 
a Peine of ne errors thus ER are SR with Be De by v. N ises ung @. 1 G h ul 2; and a Mr 
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Onenka omuÖkn, anna Tonsmunom nun EG N poTpANON Anno 
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 HWANbHEIX yYPABHeHNHÄ ePBOoro TMopsamka, Obodmaerca Mm YTOuHAeTrca ua IpmÖnnenns 
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. Eine Fehlerabschätzung für das Adamssche Verfahren Korttsränder. Et rabolattone wurde k 
‚zuerst von R. v. Mises!) angegeben. Für das Interpolationsverfahren, ‘welches durch iterative 
Verbesserung eine höhere Genauigkeit erreicht, wurde sie nach demselben Prinzip von G. Schulz?) 
geleistet. Von einem ganz anderen Gedanken ausgehend teilte später Tollmien 3) für beide 
Arten des Verfahrens und den Annäherungsgrad n—=3 eine Fehlerabschätzung mit, die sich 
sowohl in der Herleitung als auch in der Anwendung durch besondere Einfachheit und Handlich- 
keit gegenüber den Abschätzungen von v.Mises und G. Schulz auszeichnet. 
Er: In der vorliegenden Arbeit wird die Fehleruntersuchung von Tollmien auf einen beiahien. 
.: Annäherungsgrad ausgedehnt. Eine ebenfalls von Tollmien angedeutete Möglichkeit zur Ver- 
ä schärfung der Fehlerabschätzung, deren Ergebnis für die extrapolative Form des Verfahrens 
und n=3 dort angegeben ist, wird für beide Arten des Verfahrens allgemein durchgeführt. 
Anschließend wird die naheliegende Frage untersucht, welche der beiden nun vorliegenden Fehler- 


- 


Er ' abschätzungen — die von v. Mises und G. Schulz einerseits und die verallgemeinerte und 
2 .. verschärfte Tollm‘ensche andererseits — die gunsugete: Fehlerschranke liefert. Paspier 
 . erläutern die AIRERAUNG der Fehlerformeln. 

- $1. Die male Fehlerabschätzung 


u . 74% In a Ditferentialglichung a 
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ersetzt werden. 

Das Adamssche Verfahren zur numerischen Integration der Differentialgleichung (1) wird 
in der extrapolativen und interpolativen Form, wie es in den genannten Arbeiten von v. Mises 
und G. Schulz beschrieben ist, als bekannt vorausgesetzt. 


Zur Fehlerabschätzung benutzt Tollmien folgenden Satz aus der Theorie der Differential- 
gleichungen®). Gegeben seien die Differentialgleichungen 


etle); a CE RC RE er 35 


1) R.v. Mises, Zur numerischen Integration von Differentialgleichungen. Z. angew. Math. Mech. 10 

1930), 8. 81. 

2) Dr üntherSchulz, Interpolationsverfahren zur numerischen Integration gewöhnlicher Differential- 
gleichungen. Z. angew. Math. Mech. 12 (1932), 8. 44. 

-3) W. Tollmien, Über die Fehlerabschätzung beim Adamsschen Verfahren zur Integration gewöhn- 
licher Differentialgleichungen. Z. angew. Math. Mech. 18 (1938), S. 83. — Hinweise auf einzelne Formeln 
dieser Arbeiten werden kurz durch M, S, T mit der betreffenden Nummer gekennzeichnet. 

*)8.2.B. Ka Be; Differentialgleichungen reeller Funktionen. Leipzig 1930. II. $11, 52 Satz 2. 


E sei f (x, y) für den im senden in Betracht kommenden Bereich von x und y stetig und genüge _ Sl 
| ‚einer ee Be; - 
| Kay) fenlsMa—ul..::...-00:0)- “R 
4 "Setzt man außerdem voraus, daß feine partielle Ableitung nach y hat, so kann ® durch | 
; die schärfere Bedingung 

3 10 | | 

7 °y Sr tenlga: leistet ned tin nie Mahn Pl aa a Ars el Sa Höre (8) 
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BRETT a ug a Da a en 
1 np red Se TEN ; ‘E i Y TE er ws > 
N TIEREN EL aan 2 RD 
% * n ’ 4 ni u s Kr wur ve I ar: a, DS „“ ee ER 
2 RAVEN ad Liekaschredn sera Be. 2 2. 
222.9 ff y) und &(z,y) seien in einem Gebiet @ der 2, y-Ebene 
RL: überdies in @eina Lipschitz- Bedingung (2) miteinere heitl | Er 
ei dem gelte in @ die Abschätzung Se ST Ar ET En 2 AN 
Eh \ » ea 9l'se Fe EA m N = S & 3 


j Dann gilt für den Unterschied zweier Integrale der beiden Differentialgleichungen, die für 
ErFE\ gleiche Anfangswerte y,= n, haben, Fr | el. 

| h y—n]l seo ET A AR za rn. Da 
und diese Abschätzung gilt solange wie yund 7 dem Gebiet @ angehören, was nachträglich mit 
Hilfe der Abschätzung geprüft werden kann. Sind die Anfangsbedingungen der beiden Integrale 
‚nicht gleich, ist etwa Y(&,) = Yo» (0) = no und |y —no| SZ ö,, so lautet die Abschätzung?) 


vorlagen) oe SEE Eee 


Dieser Satz ist offenbar geeignet, den Fehler abzuschätzen, der bei der Integration einer Differen- 
tialgleichung y’=f(z,y) begangen wird, wenn die rechte Seite durch eine Näherungsfunktion 
fiz,y) + ®(x,y) ersetzt wird. Einige Beispiele zur Illustration dieser Abschätzung werden 
in $5 im Zusammenhang mit der verschärften Fehlerabschätzung betrachtet. 3 

Die nach dem Adamsschen Verfahren für die äquidistanten Werte «,— 2z,+»'h(h=Schritt- 
breite, v»=1,2,...) berechneten Funktionswerte y,sind bekänntlich nicht Werte der gesuchten 
Lösung y(x,) der vorgelegten Differentialgleichung y’= f(z, y), sondern sie sind Werte einer 
Näherungslösung %, die eine? angenäherten Differentialgleichung 


V=f(2.%(0) 
genügen. Um den Fehler y— %. nach (6) oder (7) abschätzen zu können, muß zunächst die 
Differenz f— f= # abgeschätzt werden. Die Funktion f (2, y) besteht aus Polynomen n. Grades. 
Bezeichnet man mit P%”(x) ein Polynom n. Grades, welches für = 2, 24.1»... Zyn die 
Werte F,=f(&, %®)) Fr --., Fn annimmt, so ist beim interpolativen Verfahren®) 
für fim Intervall Pa 8 Po x) und in den Intervallen ,_.,°,=n-+l,...) P® (2) 
zu nehmen. Beim extrapolativen Verfahren ist po x) in den Intervallen a, 2,_,»=n-+J1,...) 


zu nehmen; die Berechnung der Anfangswerte wird hier ja außerhalb des Verfahrens vor- 
genommen. 


Der Unterschied f — ei ist also der Interpolationsfehler, der beim Ersatz der Funktion 
(», y(«))=F(x) durch das Näherungspolynom P® (x) begangen wird. Dieser beträgt aber 


“ 
ep nt rn ann he nn ni De 


an) le a 4 

Ar) = GET RT RE RE F (Eh ren (8), }: 

a 

wobei & ein Mittelwert ist, der zwischen den zur Interpolation benutzten x-Werten liegt. 4 


(8) ist jetzt für die in Frage kommenden Bereiche abzuschätzen. Setzt man 


ı— n=ulh, 
so folgt 2 
ey hr+ 1 1 b; Fin+ N . 
> nEn wu +1) (u +2) .*-(u-+n) Ce 2: (9). | 
Beim extrapolativen Verfahren kommen die Intervalle z,**- Yan, n-+l,..,) 
in Frage, d.h. v=0:--+1. Nun ist aber für 0 su s1 


uurl)au +2). mn) sm tl... ... REN REHER E . 


-Ist ferner |F®+D(&)| < Ka, für den gesamten Integrationsbereich %’""%, so gilt nach (7), 
wehn mit ö die absolute Fehlerschranke der Anfangswerte 4 Ya +.» Y%n bezeichnet wird, für . 
Ur ine 


un „M En) .— hut 1 a (e" (%-%,,) CIE 1) 


ir 
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°)8. Kamke,a.a.0. I, $ 11, 51, Satz1. Der dortige Satz ist hier spezialisiert auf ,= 0, = E&. 

°) Es wird vorausgesetzt, daß die Anfangswerte beim interpolativen Verfahren in der bei G. Se hulz 
angegebenen Weise auch mit einem Polynom n, Grades berechnet. werden, so daß Y„ noch mit zu den gesondert 
zu berechnenden Anfangswerten gehört. | 


fi 


je Fch 0! e | für ji ee RE zu ı gewinnen, 
polatio spolynome Pa) ‚bei, der ‚Ermittlung der. Anfan; 
später in Bereichen Lot ip v=n+1, u ‚benutzt v 
h (9) muß daher für diese Bereiche, d.h. für ER 


8) 
1. = u<o ‚bei allen späteren Integrationsschritten abgeschätzt werde 
ipt vorkommende Bereich für u De also durch —n Ss u = 0 Sr ee. 


E Führt man in I die Substitution v=Vv Ne DR: durch, so sent 


Aa { 


er BE Wurde +2: urn 
EN im ı Falle eines geraden n=2m über in 
EEE Al) = ey en 
und im Bo ‚eines N n=2m+lüberin RN 

t a BR 
r a Die abzuschätzenden F unktionen PA )) und 92 (v ) haben n+1 einfache Nullstellen bei _-, ‚ hs a, 
E “ion + 1, “ “ 5; ; sie ‚sind ungerade bzw. ee Funktionen, Man kann sich also zur Ermittlung | i 
be Be größten absoluten Betrages auf positive v <z ” beschränken. Außerdem kann man zeigen, ER 
f Fa daß das Maximum des absoluten Betrages stetsin ir) Intervall 5 Er hd > bzw. an Er =H 1 v 


liegt, also stets ‘zwischen den beiden -äußersten Nullstellen. 


"Zum Beweis gehen wir auf u(u +1)» (u +n) zurück und setzen u+n=w. Esgenügt 
nun zu zeigen, daß die Funktion | 


hu) = w(w 1): (w—n) 


% im IntervallO---1 stets absolut größer ist als in jedem anderen Intervall v---» +1 (v al; 
7 3%... m —2). Bezeichnet © vorübergehend einen echten Bruch, so ist 
E Fe hi) Bl) (dom) zaN Pe 
N hd+r) BH+r)B+r—1) +) — 1). (d—(n —-v)) Be HR 
i a on pe 0) NW en | 2 
u: ei P+NW—1+9):.-(1+P%) PH! re 
i un 1.uwegen n—- m >22, rt; it 
= Also 
2 | | ro] w +9] 

für 


Fr RR ARR Deer RE 
RR läßt sich zeigen, daß für v <n—2 
Ä Am--1 +92] +9). 
Das gesuchte Maximum des absoluten Betrages von b (w) liegt also im Intervall 0 <w <1 
bzw.n— 1 <w <n; welches der beiden Intervalle wir betrachten, ist nach den Ausführungen 
e “über g,(v) und g,(v) gleichgültig. 
3 | Es gilt daher nach (7) mit Einschluß der Anfangswerte die Abschätzung: 


m hit 10n tı1/_M (z-s,) | 
ne“ 3 e ER LER er Di VrlLd 
i ”, Fürn=3 RN an in T (56) über. T (56) enthält einen Druckfehler: es muß ö an Stelle von y H heißen. 
Das gleiche gilt für T (56a). 
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(14) stimmt für » = 3 mit T (47) überein. 


\ 5 _ polative Verfahren wesentlich genauer ist als das,extrapolative Verfahı 
der Anfangswerte beim extrapolativen Verfahren unberücksichtigt 
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Die Abschätzungsformeln (11) und (14) zeigen die beka 


un 


Formeln nur um den Faktor Ba der folgende Werte Balz Pas css A 


Ah 1; ER NE Eh E23 Le DREH, 
0,185; 0,0642; 3 0,0417; 0,0296; 0,0232 ie 3 e 


SET: 82. Die verschärfte Fehlerabsehätzng 0... 
In der zit.' Arbeit erwähnt Tollmien eine Möglichkeit zur Verschärfung der durch 


| gegebenen Abschätzung, die darin besteht, daß an Stelle der bisherigen Definition für e 2 ar 


neue tritt: 
k 


IKIOLE: 


. Für das extrapolative Verfahren und n—=3 gibt Tollmien das Ergebnis dieser Verschärfung an. 
Diese ‘Verschärfung soll hier etwas allgemeiner untersucht und auf beide Arten des Adams- 


<E e— &| . 


% 


schen Verfahrens angewandt werden. Wir betrachten wiein $ 1 die beiden Differentialgleichungen 


y' = f(a, Y); 7 = fie, n) +9 (2, n) “ 


mit in @ stetigen Funktionen f(z,y) und 9(2,9); f(&,y)-erfülle die Lipschitz-Bedingung. Man 7 
kann die erste Differentialgleichung als die zu lösende, die zweite als die gelöste genäherte Diffe- 


rentialgleichung auffassen. Wir betrachten also die Lösung (2) als bekannt und den Fehler 
|y—n| als gesucht. Wendet man das Picardsche Iterationsverfahren auf y=f(z,y) so an, 
daß »(x2) die Ausgangslösung ist, also “ 


Ha) Yo + Ita nto))dz 


und 


Yy(X) ER Yo +/ Rs Y-1(2)) de ’ 


so kommt man in bekannter Weise zu der Abschätzung 


y—n]l sl 1) TREE INTER (18), 
in der die Konstante & durch 
z 
Jo, n(a))dx| sele— | 
% 
oder auch 
? 
OR | Plz, nl2)) BR n.. Ws ER ae (19) 
Ba a‘ 


definiert ist. 


Danach ist e nicht mehr wie bei der einfachen Fehlerabschätzung das absolute Maximum 
der Funktion 9(%,7(#)), sondern das absolute Maximum des Integralmittels dieser Funktion 
im Intervall z,°' 2. | 

Ist speziell 9(2,n7(2)) Z:0 in @,***& und monoton fallend bzw. wachsend, so liegt das 
absolute Maximum des Integralmittels bei x, bzw. x; ist 9(&,n7(2)) <0 und monoton fallend 
bzw. wachsend, so liegt das absolute Maximum umgekehrt bei & bzw. z,. Ist Yla,nla)) in a: *® 
von wechselndem Vorzeichen und monoton, so liegt das absolute Maximum des Integralmittels 
entweder bei x, oder bei x. Ist schließlich 9 (&,»7(&)) nicht monoton, so liegt das absolute Maxi- 
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mum im Innern des Intervalles ©, 2, und zwar an einer Stelle"&, für die der Funktionswert 
gleich dem Integral ist, denn aus u 


folgt 


u ER 
5 % een | PLRPER ER  E EE N En 


n 


Bei der Anwendung auf das Adanhssche Integrationsverfahren tritt noch die Vereinfachung 
ein, daß ® nur von & und nicht von y bzw. n abhängt. 

Die Anwendung der verschärften Fehlerabschätzung wird durch folgende Eigenschaft des 
Integralmittels erleichtert. Wir betrachten eine nirgend negative periodische Funktion f(z) 
(Periodenlänge A), deren einzige Nullstellen beivA #=0, +1, 42, ...) liegen, außerdem sei f(z) 
im Innern jedes Intervalles »4---(»+1)4 differenzierbar und habe dort genau ein Maximum. 
Das Integralmittel dieser Funktion ER 


Ja) — 2 [fta)dz 
0 


hat nach (20) Minima und Maxima an Stellen £&,, die sich aus der Bedingung ergeben 


3 & B 
1 
eg BR N RR N a rt 21). 
fi) 5] fta)de | 1) 


Dort also, wo die Kurven f(x) und J(x) sich schneiden, liegt ein Maximum oder ein Minimum 
von J(z). Da 
FE) 


&, 
haben f’(x) und J’ (x) für &,>0 gleiches Vorzeichen. Dort also, wo J(x) den aufsteigenden Ast 
von f(#) schneidet, liegt immer ein Minimum, dort, wo es den absteigenden Ast schneidet, liegt 
ein Maximum. 

Weiter folgt aus (21), daß die Maxima mit wachsendem & beständig fallen und die Minima 
beständig wachsen, : 


I) = 


TR 


2 
denn das erste Maximum liegt auf dem absteigenden Kurvenast bei &, und ist durch die 
Gleichheit der Flächeninhalte « und 5 bestimmt. Das zweite Maximum kann nun aber nicht 
die gleiche Größe haben, da dann das Rechteck 0, (1), (2), (3) (s. Bild 1) einen um c größeren 
Flächeninhalt hätte als f(x) bis zur Stelle (1). Also muß das zweite Maximum kleiner sein als 
das erste und so fort. In.ähnlicher Weise kann das Anwachsen der Minima klargemacht werden. 
Schließlich gilt, wenn n > 0 ganz und O<d<I1 


N 
7 


Bild 1 Bild 2 


(n+ 8). \ mat dA 
1 
ums) = im 5 j fi2)dxz= lim ler + j fd) 
To von (n a d)2 Pr No (n +D)A nA 
N a 
ag JA): mit PF=/[f(a)de, 
v 


so daß J(x) eine um die Gerade y = Z mit abnehmenden Amplituden pendelnde Kurve ist. In 


A 
dem Fall, daß f(x) in jedem Intervall monoton wächst, sind alle Maxima gleich F/A und liegen 
an den Intervallenden (s. Bild 2). 


| ne ee . ’ 
M=(s— 1%) ee 
gilt. Weiter ist, ha die Re bei A er ae Lu 
se 


wenn 


| / De el SKaı | 
für den gesamten Integrationsbereich %o’ "© gilt. Nün ist a in a ne Intervall, u Sa 


monoton wachsend, so daß der in Bild 2 dargestellte Fall v demnach: 
® a4 r Pf) 
1 . u! 
= Pads — | u-me-n) ae ne| , 0 
2 — I Ant Am) 7 A 
| tt [ulm +1)--- (u +n)du=h"t!(n RN Ä % | 
wenn man 2.— m = u und die in der Integrationsformel für das extrapolative Verfahren 
auftretenden Koeffizienten Ay £ 
nun +1)-- Er FE et, 5 (24), 
die bekanntlich sämtlich Ne sind Se = 1/2, = 5/12, a,— 3/8, .. .), einführt. Also ergibt (22) 
Ray u jpau SR (a PORT R (25), 
und die verschärfte Fehlerabschatzung für das extrapolative Verfahren lautet nach (7) und (19) 
v—9l SUEN, EURER ui er ha N EN: 


Der durch die Anfangswerte bedingte Fehler ö bleibt von der Verschärfung unberührt. 
Der zweite Summand in (26) ist also verglichen mit der einfachen Abschätzung (11) um die 
Koeffizienten a,<1 verbessert worden. Für n=3 geht (26) inT (56a) über?). 

‚Führt man für das interpolative Verfahren dieselbe Überlegung durch, so ist wegen (8) 


Sa 3 )de Ss — IE ger Firolas DENE (27) N 

mit 2 
NR ee Y-ı Ssesy Pp=en+lınr2,...) Is) 

m = (0 a) (2): (20) u Su <a 
Aus den Überlegungen von $1 folgte nun, daß (x) im Intervall z,*** x, und %,., "2, an i 


jeder Stelle größer ist, als an entsprechenden ‚Stellen der anderen mar: zwischen x, und 2,. 
Außerdem ist der Verlauf von m (2) in den Intervallen ab x» derselbe wie im Intervall Sa Ya 
F ührt man daher die Funktion QP’(«) ein, die in allen Intervallen mit Ausnahme des ersten 


u u EA nn nn nn 


ee wege 


%,"*"%, denselben Verlauf haben soll wie la” w)| in au. *** Zu bzw. den späteren Inter- _ E 
vallen, während sie im ersten Intervall mit |q%’(«)| übereinstimmen soll, so gilt Ku 3 
| ga) | So (a (2) $ 
und 4 
Pt | Ä 
IN a | WOLLEN ER | 

a N ih: = —%(r +1) (n+1)! On Ey) 
y ..’® 
®) Vgl. Fußnote 7, !l 


Be rin 2 + 


itervall erreicht, denn wenn 


er 
r . 


an A, 
m > y Y Ne ne i i 


‚ab, das ereie, Maximum also Bei & Hegt,.s so Kaun dag Br ol die gleiche Größe rs 


BE 


ar so daß gilt u e de | EA 2 16,910; 
Io a = . jeo dx wen Jene a ... a asıaıa 
i WERD 
6% Fe < Max fu & ja (u ... ra mia) ee en Bi: 
wenn 2 — = uh gesetzt wird. Demnach geht (29) über in : 
Labs „ Jroaz|= Set Ani ax r ee nt (u—n)|du R 130) N 
7% (n we Di ' Rs = 
E ! LE ‚ Dr 
Setzt man noch PN U Mc Br | 
{ # nr. u“) N Be a P i 
ve , Nus(t j A uonldu) - SEE BEN NE 
so lautet die verschärfte ee für das interpolative Verfahren er $ ä 
13 AT Kr+i M (0-%) —_ 1 i | | 
v9 Si u‘ N ru ti 
In der Fehlerabschätzung (14) stand e„ an Stelle von ©,. Für (15) kann auch 
&n = Max |Ju (u— 1) (u—n)] 0<u<sl 


geschrieben werden, da 28 Maximum von u (w+1) (u +n) für -1su<s0 dasselbe ist 
wie das von |u (w— 1): -(u—n)| in su s 1. Jetzt tritt an Stelle von e„ das Maximum 


des Integralmittels O, < en, so daß - — <1 die Verbesserung in der Fehlerabschätzung bedeutet. 
Die folgende ZaRele enthält eine een bartalte dieser Werte: 


AS R 
0,385 1 FRE 55 
O, 0,188° 0,296 0,776 SE Re} 
= 0,75 0,768° 0,776 0,800 


Aus den gegebenen Fehlerabschätzungen (11), (14) bzw. (26)» (32) folgt, daß der Fehler von der 
(rn +1)-ten Ordnung in A ist und daher 
j h->0 Er 


so daß die Konvergenz des Verfahrens gegen die wahre Lösung gesichert ist. 


da 2c>d,.denn es ist c vi 6 Dee .d. Die folgenden, Pass Patien wieder eine > fallende e 


Die o @,,, sind die in der BE ERHSN.SN für das extrapolative Verfahren 


2 ist die u reelle Wurzel der Gleichung 


Die von v. Mi Mises gegebene Fe 18 r das e ii olatıve \ rfähren 
De hier Be ne GC) ET E. > ER u 
| ur: 
N x . u _— u Bat ). 
” W—7] <or Hm ee -n: En = ER ; 
ji Darin ist » die Schrittzahl, VER, | ERLINR: 


“ 


N {a R* Redman 
zienten, sofern man nach den SEN RARE NER F; an Stelle der es abe er 


sie gilt | FR Ir. Y 
. '@,, ey D3 ER je * .. ae NE iR 


W 


r 


Pe D)=EME ana AR ee N ee | 


die >1 ist. Die 2-Werte können für eine Reihe von n- und A M-Werten aus den ERENT 0; 
stellungen bei v.Mises näherungsweise entnommen werden. IR 


Es sollen. die Abschätzungen (34) und (26) miteinander verglichen werden. Man muß dazu 
bedenken, daß 
C— I k \ 
h 


Ye 


die Gesamtzahl der Integrationsschritte und 
C— In 
en ar 
die Zahl der Integrationsschritte von x„ ab bedeutet, so daß n +u =» gilt. 

Das Verhältnis der beiden ersten Summanden in (26) und (34) ist dann 


und es ist für alle u IE 1% sofern z>ehM ist. Für h—0 folgt v=1, da z=1. Durch Differen- 
| tiation von (37) unter Be- 
rücksichtigung von (35) er- 
gibt sich 
dz 
ur TAM Au - (39). 
Da aber, wie aus (36) folgt 


B> ı=1 . (40) 
i=0 


und daher A„>1 sein muß, | 
haben die Kurven z=z(kAM) e } 
für AM=0 einen größeren . r 
Anstieg als z=ehM, so daß 

für hinreichend kleine 


Bild a hM z>e*Mist. Da anderer- j 
seits, wie aus (37) ersichtlich, h 
h M mit wachsendem z.von erster Ordnung unendlich wird, MUB für hinreichend große RM z <etM ’ 
sein. Die numerische Berechnung ergibt ; 
$ ulNn; ar 3; 4 | 
hM » ’ 
e>e Für AI ARTE 1,46; 1,78; 20ER (Al). 


Da die praktisch vorkommenden h M-Werte weit unterhalb der Schranken (41) liegen, ist tat- | 
sächlich stets » > 1, und zwar um so größer je größer 4 wird. 

Auf Bild 4 ist der Quotient p# für „= 1 und u=20 dargestellt. Man sieht, daß der 1, Sum- 
mand der Abschätzung (34) um so mehr das entsprechende Glied der Abschätzung (26) überwiegt, 


TE N NE 


/ 
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je größer AM, nund usind. Insbesondere wird die von Misessche Abschätzung für große u 
praktisch nicht mehr brauchbar, da die Fehler um Zehnerpotenzen größer werden können als 
nach (26). Hinzu kommt, daß der Quotient 9" noch mit z">1 zu multiplizieren ist. 

Für das Verhältnis der 2. Summanden von (26) und (34) gilt | 


A zutnr — 1 \ 
Ae = An Mm] N USE A Te Ra (42). 
Dann ist, wenn (41) vorausgesetzt wird, Wi 
Bm Aa a REN ee SR (43), 
HU>o 


d.h. für große Schrittzahlen u wird auch das Verhältnis A, beliebig klein. 
Für A—0 ergibt sich 


REST NER ea 


h—0 va—o dz 


und wegen (39) 
; al 
se 
h—0 n+tu . 


Die numerische Berechnung ergibt, daß auch für alle A>0, A,<1 ist, sofern die Bedingung (41) 
erfüllt ist. | 
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Bild 7 Bild 8 


Auf den Bildern 5—8 ist das Verhältnis A, für n=1,2,3,4, u s30, kM <0,2 dar- 
gestellt. Die Kurven zeigen das (43) entsprechende Abklingen für große u. Sie liegen für alle 
betrachteten n und AM unter der AM= 0 entsprechenden Kurve. Das Verhältnis A, wird außer 
mit wachsender Schrittzahl auch mit wachsendem Grad n und wachsender Schrittbreite immer 
günstiger für (26). Wegen » > n sind diese Kurven erst ab «= 1 zu betrachten. 


EA FE NEE 
‚ WEREUESK, + Ben 
“ Pr Mar 


der 2, Summandın (Hit 


Glieder von (6). fi N n pr: 
Insgesamt ist also die Abschätzung (26) ‚der isch immer erfüllten Be 
für alle m und w günstiger als die Abschät (34); insbesondere wird (26) für-gı 


. zahlen von geringerer Ordnung unendlich als (34). Für nach Null gehende Schrittbreii 


beide Abschätzungen für große Schrittzahlen ineinander über. hl Tee En : 


$4. Vergleich der Abschätzungen für das interpolative Verfahren 

Nach G. Schulz (43) gilt für das interpolative Verfahren TR 

| NSS Butıl Entı ,_ Rn. 
| y—y| ae vr Bel Tu er nr es er: 
RE DR f r>n a u 


e e N 
Ar= 5 |bn,.| =. win ae TE 
: . 4=0 i ; 3 . :- 
Die bn,, sind die den @,2 entsprechenden Koeffizienten des interpolativen Verfahrens, für de 


mit 


n # 
 : tab) 
a \ x ' 
gilt, (40) entsprechend folgt hieraus i ö 
| De DE RR (47), ° 
b i=0 
so daß auch AR >1 für n>1 ist. z ist die kleinste reelle Wurzel der Gleichung 
PAR —DEIM End)... a Be (48), 
1=0 


die’>1 ist. | | 
v ist wieder. die Gesamtschrittzahl und e eine Schranke für die absoluten Beträge der 
Fehler &,&, -=., &m, der Werte y(z,), Ya) --, Y(za) 2). 


Obere Schranken für &,, &9, . . .,&n hat.man nach S(72) in den positiven Zahlen Y,,, die aus 
den Gleichungen 


$ f 
Eu ala Bei u j 


zu bestimmen sind. Da Y,>d,, soll zum Vergleich der Fehlerabschätzung (32) mit (44) an 
Stelle der Y, d, genommen werden. Die Größen d, sind definiert durch S(74): 


k 
v | 
y=hnt2 Kur? 6 =) AN RR N ‚ (49). | 
k-llk—ı I 
Führt man die in der Integrationsformel zur Bestimmung der Anfangswerte auftretenden 3 
Koeffizienten i; 
{i) f 1 
1 A 
Bm = 77 (“—m) u—m +1) u—m+i—1)du ...... (50) » 4 
| | En | 
‚ein, so geht (49) über in 1 

y 

Sa Raı & |Br-ı,n+ıl She: ER I (S1). i 
Es kann demnach für e der Wert d„ genommen werdön. Das Verhältnis 4; der Fehlerabschätzung 1 
(32) zu (44) ist dann . 1} 
vyhM | 
HB ee 50 H 
‘ &nhM2+Pß, (2? E57 1) BR ee AR ( ) i 
| 


en # 
®) Da hier (vgl. Fußnote 6) y,„ mit zu den Anfangswerten gehört, ist & als Schranke der Werte Eur Een 


und nicht wie bei G.Schulz für &, &,...,&,_—. aufzufassen. & und &, sind von den ebenso bezei 
Größen in $l und $2 zu unterscheiden. 0 > nn SENEr 


un TE Re 


Y 
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mit 5 Ei Be 
n+N!X n+l! nzıl 
= — B I Re ——— TE... nn 93): 
n O, 2 | 1,1] Pn @3 002048 (93) 


‘Man findet: 
n | 1 2 3 BES 
&n | 0,885 1,69 2,10 DO ler (54). 
BR’ BI? EO:BB5, 150,798: P010,5B7U, 0442, | 
a Auch hier ist 


v—&© 


TE ae ee a (55). 


- einen Pol haben, muß 
n,Nn 


’ 


Da die Kurven z—=z(hM), wie aus (48) folgt, für AM = 


Da außerdem entsprechend (39) 


ESEL NEN ER ER EN. Ye a u) 


> dM 


Ba 

Fa NR B2, 

für kleine AM die Bedingung z>e%M erfüllt. Für n=1 haben beide Kurven den gleichen 
Anstieg. Eine Reihenentwicklung der a=1 entsprechenden z-Kurve 


_gilt, ist auch wegen 


2 2 M)® hM): 
un m ‚oa ame. 
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70 202 y— 3% 
Bild 10 


N 


Bild 11 Bild 12 


n q 
für große AM z>e*M erfüllt sein; wegen der Konvergenzbedingung S(13) bleibt AM <1/a,. 


ist, ut RR aut i RR Ss müßt 


ein mindestens zweimaliges Schneiden der z-Kurven mit e®M# 
Im Gegensatz zum ne Verfahren ist also 2> ehM 


'  Schrittzahl wird die Fehlerabschätzung (32) in 


in jedem Falle günstiger als (44). 
Auf den Bildern 9—12 ist das Verhältnis 4; für n=1,2,3,4 hM 0 v <i 
wachsender Schrittzahl », 


gestellt.” Man erkennt das Absinken von A; mit er 
ne so eher erfolgt, je größer hM und » sind. Verglichen mit den.Bildern 5-8 sicht man, er 2 


die Unterschiede beider Fehlerschätzungen jetzt geringer sind. . 
Für verschwindenides hM erhält man wegen ch b 


yeei FR = 
Se ade 77 ER 2. 


Auch dieser Grenzfall ist auf den Bildern 9—12 dargestellt. Da alle Kurven unterhalb 


dieser Grenzkurve liegen und da diese Grenzkurve für »—>oo dem Wert BT ki: 


1 1 
Bade (nr +1)! BR 
zustrebt, kann ER Verhältnis A; für sehr kleine AM vorübergehend etwas größer als 1 werden; 
wie auch aus den Bildern 9—12 ersichtlich ist. Wegen » > n sind die Kurven bis »—= n gestrichelt. 
Ein Zahlenbeispiel: Für n=4, AM=0,1, v»=50 nach Bild 12 4,=0,07, d.h. Abschätzung 


(32) ist rund 14mal besser als (44). 
Die verallgemeinerte Tollmiensche Fehlerabschätzung, die sowohl in der Herleitung als 


‚auch in der Handhabung einfacher ist als die entsprechende Abschätzung von v.Mises und 


G. Schulz, erweist sich für das extrapolative Verfahren unter der praktisch immer erfüllten 
Bedingung (41) stets als wesentlich günstiger und für das interpolative Verfahren mit Ausnahme 


. sehr kleiner k M-Werte und beschränkte Schrittzahlen, für die beide Abschätzungen die 


Größenordnung besitzen, als günstiger als die Fehlerschranke von v.Mises und G. Schulz. 
Insbesondere wird für das extrapolative Verfahren unter der Bedingung (41), für das- inter- 
polative Verfahren stets für große Schrittzahlen der Unterschied beider Fehlerabschätzungen 
beliebig groß. 


$ 5. Beispiele zur N. Fehlerabsehätzung 
Um die Güte der gegebenen Fehlerschranken zu prüfen, soll zunächst die grundlegende 
Beziehung (6) mit (5) und mit der Verschärfung (19) an drei Beispielen behandelt und mit dem 
wahren Fehler verglichen werden. Im Hinblick auf das Adamssche Verfahren soll 9 nur von & 
abhängen. 
1. Beispiel: Es seien die Differentialgleichungen gegeben: 


RR Lars r BR 
TEEN TS EN (58). 
Die gemeinsame Anfangsbedingung sei 2,—= y, = 1. Gesucht sei der Unterschied |y— | bei 
x=2. Als Gebiet G kann der Parallelstreifen 1S2<2 betrachtet werden. Man kann M=1° 
und nach der einfachen Abschätzung, da 9(2) = « ist, 
da)| S2=e 


nehmen. Dann gilt nach (5), (6) 
gen ste —n), 
Nach der verschärften Abschätzung (18), (19) ist 


© 
EZ - Max (+1) 2, 
1, 2 2 


y— nl <a (@1—1). 


1 
x 
so daß 
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Der wahre Fehler ist SER, | ; | Br 


ynl=#—e, 
so daß der Fehler f bei »=2 in der Reihenfolge einfache Abschätzung, verschärfte Abschätzung, | 
wahrer Fehler folgende Werte Annimmt: 
EAN E258: vn 


2. Beispiel: Während im 1. Beispiel 9(x) monoton wachsend war und das Maximum des 
Integralmittels demzufolge am Ende des Integrationsbereiches lag, soll jetzt ein Beispiel mit 
nicht monotonem #(x) betrachtet werden. Es sei 


Ve BEN en RR TE Be (39). 
Anfangsbedingung sei &,=0, %—=1. Es kann wieder M=1 genommen a Gesucht sei 
der Fehler bei z=n. 
Als Gebiet @ kann der Parallelstreifen 0 <x <n betrachtet werden. Da 
O@)I= Isin.z| <1, 
gilt nach (5), (6) 
y—alse—1. 
Nach (19) ist \ 


ee j sinzdz )- Max ( — — 0,725 
x & 


0 


a 


also 
y—n| s 0,725 (ee —1). 
Der wahre Fehler re sich zu 


y— 1-— ze + 6inz +FC0S'%), 


so daß der Fehler bei = sich zu 
RN De ER 
ergibt. 
3. Beispiel: 
ENDEN TEN TEEN N ra (60) 


mit der Anfangsbedingung y(0) = 1. Wir betrachten das Intervall 0 <x s1. In diesem Be- 
reich ist 9—=1— x monoton fallend, das Maximum des Integralmittels also gleich dem Funktions- 
wert für ©=0, also e=1. Da auch M=1, gilt nach (18), (19) 


Bunzse—!: | 


Die einfache Fehlerabschätzung fällt natürlich mit der verschärften zusammen. Demgegenüber 
ergibt sich der wahre Fehler uf=z, day=e, n=x-+e. Also gilt für =1 


fsıR; fl. 


$ 6. Ein Beispiel zur Fehlerabschätzung für das Adamssche Verfahren 
Als Beispiel sei wie bei v. Mises und G. Schulz die Differentialgleichung 


MENT EURE Den Kerr (61) 


mit dem Anfangswert 2=0, y=1 mit h=0,02, n= 4 gewählt. Für das extrapolative Verfahren 
seien die Anfangswerte in der dort beschriebenen Art gewonnen, so daß ö= 0,5 - 10”® gesetzt 
werden kann. Nach S.p.52 ergibt sich mit (34) nach v»=100 Schritten die Fehlerschranke 
l&ıoo | < 1990 - 10°®. 
Nach der Abschätzung (26) erhält man dagegen mit (s. M. p. 90—91) 
hM=0,006; M(z=— x,) = uhM = % - 0,006 = 0,576; 
95 


K, = 8840; 4 = 988 == 0,3299 
8840 0,576 £ 
leioot S 0,5 : 107° - e0:578 +.0,3299 - 0.006 « 64 » 10-12 (60,576 — 1) = 25 - 10-8 


also eine rund 80mal bessere Schranke. 
12 


PER, 


ER 


4 


Br er 2 ee 


N so wiel Stellen dı geführt is 
MR ‚die Anfangswerte auch mit einem eh rad 


Ko SER E hat n ZBAR Da en 


RR iR: er 
RE, i d - 28. 10-8 
Eine einmalige Iteration mit Hilfe der Gleichungen“ (8 ey en 
ER 0 94=3,02-10-*, ar) 


MR 


g Nimmt man diesen Wert für e und die anderen Werte wie in der Arbeit von G. Schulz, so 


für das interpolative Verfahren. nach En ERTL N Yes 


I“ 


Die Abschätzung (32) liefert mit denselben Werten wie oben: 


Kader 64 - nn . 8840 s, 
also eine — wegen yo kleinen Se „Wertes — nur ae biiwere REBEL! 


' Die Frage, wieweit die Fehlerschranken von dem tatsächlichen Fehler abweichen, 1äßt u 


% 


° 
nicht allgemein entscheiden, da es wesentlich davon abhängt, wie stark 2£ un und f*+» (2 vo) 


sich im Integrationsbereich ändern. 
Eingegangen: 23. 9. 1948. E 


Ben älune der allgemeinen, geschlossenen sowie auch 
| offenen Kegelschale durch Belastung ihrer Spitze") 
Von @. Sonntag in Vohburg (Donau) 


Es werden Bedingungsgleichungen zur Ermittlung der Spannungen in einer allgemeinen ange 


schlaffen Kegelschale abgeleitet für eine Belastung ihrer Spitze durch Kräfte und Momente. 


\ 
The conditional equations for Ihe determination of the stresses in general thin cone slabs while loaded 
on the top with forces and momenis are derivated. 


Les &quations conditionales pour determiner les tensions dans des minces cosses cöniques, en general, 
&lant chargees aux cimes par des forces et des moments sont d£rivees. 


Buisoparea ypasnenus 19 OnpegeneHug HANpSBKeHHÄ B Öe3MOMEHTHON KoHycHoli 06010uRe, 
Harpy:KeHHON B ee BepImmHe CHIaMM MH MOMEHTAMN. 
I. Ermittlung der Bedingungsgleichungen 


$1. Die Beziehung zwischen Spannungsfunktion und Spannung 
Wir legen auf den Kegelmantel ein Koordinatensystem s,y mit seinem Ursprung s= 0 


in die Kegelspitze. Jeder Punkt ist dann festgelegt durch seinen Abstand s von der Kegelspitze 


und dem Winkel y, den der Strahl s mit einer beliebig als y=0 festgelegten Erzeugenden des 
in’eine Ebene abgewickelten Kegelmantels einschließt. Da der Kegelmantel als frei von äußeren 
Lasten angenommen wird, und nach Voraussetzung von der dünnen Schale keine Biegespan- 
rg werden können, muß bekanntlich die Tangentialspannung o, überall ver- 
schwinden 


a a N 


Auf den einfachen Nachweis mit Hilfe der Gleichgewichtsbedingungen an einem durch zwei 
benachbarte Mantellinien und dazu orthogonalen Tre herausgeschnittenen Schalen- 
element soll hier verzichtet werden, 


‘) Für die Anregung und Förderung dieser Arbeit möchte ich meinem hochverehrten Lehrer, Herrn Prof. 
Dr. phil. L. Föppl und Herrn Prof. Dr. B. Baule meinen verbindlichsten Dank aussprechen. 
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% ne wi VERNETE NUN B Er“ 


RER 


INENSSONER TTS TPUST PORN INUDTERR EN 
{ 


"öppl [1]2) nachgewiesen hat und an Hanc Eching: leicht m 


n SH lassen sich die Be 0 und die ae Ten Uuteh? 


| ue ar, 
umgeht ern 


a worin fı und ER ren in y allein bedeuten. 6 bezeichnet die gleichmäßige. > 


Wandstärke der dünnen Schale und c,, C, sind Konstanten. 


, Die Funktion fı N der Belastung der Kegelspitze durch eine Kraft, und derreine Längs- 27, 


iptzmstind nimmt proportional mit sab. Die von der Funktion f, bestimmten Spannungen 


' sind umgekehrt proportional dem Quadrat des Abstandes s von der Kegelspitze, die durch ein 
Moment belastet wird [1]. Die Aufgabe, für eine gegebene Belastung der Kegelspitze die Span- 
nungen zu ermitteln, isti im allgemeinen statisch unbestimmt. Ihre Lösung soll im foigenden 

aufgezeigt werden. 65 


$2. Ermittlung der N f für die durch Kräfte an ihrer ir 
belastete Kegelschale 


Schneidet man. die Kegelschale längs einer Kurve s= const. auf, so muß die aus a Span 
nungen resultierende Kraft das Gleichgewicht zur Last ® herstellen. 


Um diese Bedingungen analytisch zu erfassen, denke man sich in die Kegelspitze ein recht- s + ! 


winkeliges x, y,2-Koordinatensystem gelegt, das zunächst ganz willkürlich gerichtet sei. Be- 
zeichnet man mit Z,,C,,, die von einer Mantellinie des Kegels und der x-, y- und z-Achse 
eingeschlossenen Winkel, so fordert das Gleichgewicht: 


B=JPEHPI HR: 
; 2% ; : 
mit P,=söfo,-cosl,.dy, P, Se ed dp P; [0 Er E @): 


y=0 y=0 

Darin ist 2 der Keilwinkel des in eine Ebene abgewickelten Kegelmantels. Die Winkel £ Ey 
und Z, sind Funktionen von y, die nach gewähltem x-, y-, z-Koordinatensystem bei einer ge- 
gebenen Kegelschale eindeutig bestimmt sind. Für die weiteren Überlegungen genügt es, die 
auf die Achsenrichtung entfallenden Komponenten P,, P,, P, der Last ® für sich zu betrachten, 
da. der Gesamtspannungszustand durch Überlagerung der" .aus diesen Lastkomponenten sich 
ergebenden Spannungen erhalten werden kann. Es werde die Untersuchung für den jetzt als 
Einzelkraft betrachteten Lastenteil P„fortgesetzt. Ist o,(,), der Anteil der Spannung o,,. der von 
der Last P, hervorgerufen wird, so lautet die Gleichgewichtsbedingung 


es g) rw a a EP ee re (4a) 
und. nach Einsetzen von o,(,, nach GI. (2), wobei f,—= 0 zu setzen ist, folgt 
Pose; Fr VERS LE a ER ES (4b). 


Die Integration ist hier wie in Zukunft immer über 9%. «zu erstrecken. Weiterhin darf die aus den 
Spannungen o,(., resultierende Kraft keine Komponente senkrecht zur Last P, aufweisen, was 
nach Gl. (3) durch die Bedingung 


IKA®) ° COS &,dy = 0, IKEA) - cos L,dy = VE (5) 

ausgedrückt wird. Drückt man in Gl. (2) c, durch die Gleichgewichtsbedingung (4b) aus, so folgt: 
P, Fi 

= — — 1. ——....... 0: ER ON 

I; (2) sö Eh cos, dy (6) 


Zur Vereinfachung lassen wir in Zukunft den Index an’ f fortfallen. Von den unendlich vielen 
Funktionen f und damit verknüpften Spannungszuständen, die die Stetigkeitsbedingung und die 
Gl. (4) und (5) befriedigen, wird sich derjenige Spannungsverlauf einstellen, bei dem die Form- 


änderungsarbeit zu einem Minimum wird. In den Ausdruck für die auf die Raumeinheit bezogene 


Formänderungsarbeit [2] 


2 
VE) 


wird die Längsspannung nach Gl. (6) eingesetzt. Nach Zusammenfassung aller von y unab- 
hängigen Größen in den Beiwert c ist die Formänderungsarbeit in einem von den Schnitten 


2) Die Zahlen in eckigen Klammern beziehen sich auf das Schrifttum am Ende der Arbeit. 
| - 12* 
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w 


au» Aha TE N 3 ur „HM: 
TR BR ad a EA a ed as, 
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180 Sonntag, Beanspruchung der geschlossenen 


Sr const. und s-+ds= const. begrenzten Streifen durch 


j° dy, Pe (8) 

m ee a el 5) Sat a ae VE RE ee 

AT TTF eos Gay? | 

bestimmt. | A \ 
Nach den Regeln der Variationsrechnung ändert man die Spannungsfunktion in seinem 
Verlauf ein wenig ab, indem eine im Intervall von y=0 bis y=2a ganz beliebige, aber stetige 
und endliche mit dem sehr kleinen Faktor e multiplizierte Funktion »(y) zu f h ddiert wird. 
Bekanntlich muß für die gesuchte Funktion f die partielle Differentiation der variierten 
Formänderungsarbeit nach e für e=0 verschwinden [3, S. 594], und es folgt, nach einfacher 


Rechnung: | 
[F eos ady- [In dp — [S?dy: [meos&,dy— 0 ER 
Gl. (9) muß für jede beliebige Wahl von n) verschwinden; diese Bedingung ist aber nur für 
feel Bi ee aa 
erfüllt. 


Gl. (10) in (6) eingesetzt gibt den Spannungsverlauf o,,) für jede beliebige Kegelschale 
als Funktion des Winkels £, an, den die Erzeugende der Kegelschale mit der äußeren Last P, 
einschließt: 


Re BR 
sö J cos? C, dy ' Er 


Die Rechnung wurde unter der zunächst stillschweigend vorausgesetzten Erfüllung der Gl. (5) 
durchgeführt. i 
Nach Einsetzen von Gl. (11) in Gl. (5) folgen als weitere Bedingungen die ersten beiden 
Gleichungen von (12). h 
[eos &zcost,dy=0, fecost,cost,dy=0, [cost,eost,dy=0 .. (12). 
Die dritte der Gl. (12) wurde der Vollständigkeit halber schon mit angeschrieben. Sie erhellen, 
daß die Richtung des in die Kegelspitze gelegten rechtwinkeligen Koordinatensystems nieht will- 
kürlich gewählt werden darf. Durch eine einfache Betrachtung, die wir uns hier wohl sparen 
können, läßt sich nachweisen, daß bei Vorhandensein von Symmetrieebenen die gesuchten Ko- 
ordinatenachsen senkrecht auf diesen stehen müssen, um GI]. (12) zu erfüllen. Da die Lage des 
Achsenkreuzes durch zwei Richtungen bereits eindeutig bestimmt ist, genügt es immer zwei der 
Gl. (12) zur Berechnung hinzuzuziehen, die dritte Gleichung ist dann auch erfüllt, wovon man 
sich zur Kontrolle leicht überzeugen kann. : H 
Danach ist die Last in diese Richtungen in P,, P,und P, zu zerlegen. Die Spannungsanteile 
0, (2)5 Osly) Is) Werden aus Gl, (11) nach entsprechender Vertauschung der Indices gefunden, 
wonach sich der endgültige Spannungsverlauf durch Superposition 
ergibt. Dieser Lösungsgang kann auch ohne Einschränkung auf die 
offene Kegelschale angewandt werden, da bei einem Schnitt längs 
einer Erzeugenden keine inneren Kräfte frei werden. 
Besondere Bedeutung kommt dieser Überlegung einer Kegel- 
schale mit einer seitlichen Öffnung im Kegelmantel zu (Bild 1). Da 
für den betrachteten Belastungsfall in der biegungsschlaffen Schale 
weder Tangentialspannungen o, noch Schubspannungen o,, über- 
tragen werden können, müssen die Erzeugenden des Kegelmantels, 
die die Öffnung überstreichen, spannungslos sein. Die Schale ist dann 
als offene Kegelschale, d.h. unter Ausschluß des die Öffnung eingren- 
Bild zenden Sektors zu berechnen. Zeigt sich dabei, daß die Spannung o, 
ER antnnitgener ,. dlepı die Öffnung tangierenden Erzeugenden verschwindet, so ist die 
Lösung als streng anzusehen. Andernfalls würde sich beim Über- 
gang zur ersten Erzeugenden, die die Öffnung überstreichen, eine Unstetigkeit ergeben, die 
wegen der Formänderungsbedingungen nicht möglich ist. Es werden dann unter Auftreten von 
Biegung auch in den nach Voraussetzung spannungsfreien Mantellinien Spannungen auftreten, 
die jedoch rasch abklingen werden. 


a 


I;() = 


.$4. Ermittlung der Spannungsfunktionfüreine Kegelschale mit Momentbelastung 


an der Spitze - 
Schneiden wir den Kegelmantel in einer Kurve s=const. auf und legen den Momentenpunkt 
in die Kegelspitze, so ist das Momient der Längsspannungen gleich null,und das Moment der 
Schubspannungen muß Gleichgewicht halten mit dem an der Spitze eingeleiteten Moment. Das 


F worin UBS Ken Hebel der Kraft 2 “ BR Er in Be auf die -Achse DAR Die Inte- y 
 . gration ist hier wie in Zukunft über den Keilwinkel 2& der abgewickelten Kegelschale zu er- 
strecken. Damit die ‚Spannungen Tuy@) Tayly) Teyin des in die z,Y, 2-Achsen zerlegt en 
Momentenvektors ‚superponiert werden können, muß das Achsenkreuz so gerichtet sein, daß 
BR Re Schubspannungen Toy) kein resultierendes FERN in. bezug auf die Y Reh ; 
aufweisen, woraus folgt: ER 


Ba 68? een Bd, En? a: 4% ey, [tar 0, dy=0 Pt; 1. 


"a e Für M, und M, sind die Indices entsprechend zu vertauschen. Die y 2) worin jetzt ar 0 zu R, , 
Er, ‚setzen ist, lauten in Verbindung mit Gl. FAR » Ki 


Er H | # Lee BA En x E Tl ni | ” 4 
Br a . v 7 « u, ee ee Yo Zn Y . €) 15 ’ 


worin Striche an h partielle Differentiation nach % es Die Bheiche Variation der Form- ä 
Sg en wie in $2 führt auf die as 


4 (®: aufn: dy—[fr-dy[Bn- dy) + JS Bayfn- v—SFayJOndv)=0 (16a), 


wobei auf die Bi an f und ® verzichtet a nist die beliebige Variationsfunktion BE E,6G. 
der Elastizitäts- bzw. Schubmodul. ‚Damit n nicht in seiner Ableitung n’ auftritt, integrieren 


E: . Wr Bel: N \ 
3 | If mr dy = Une—[fn- w) 
und außerdem zweckmäßig Jf?dy =[f'fE*— [ff -dy | 
Bei einer geschlossenen Kegelschale, d.h. jede Kurve s= const. läuft in sich zurück, verschwinden ° 


die eckigen Klammern in Gl. (17) wegen der Stetigkeitsbedingung, und es folgt für die erste 
Klammer in Gl. (16 a 


a a a an an 


P\:.; i = (— [r® ap [f"n - dp HF dp Dn-dy) ee (16b). 
j _ Für beliebige zulässige Funktionen 7 ist Gl. (16) erfüllt, wenn BR 
=D 01:5 = const, "idder Belt HN anne - (18) 


f . gesetzt wird, wie man sich durch Einsetzen in Gl. (16) sofort überzeugt. ® = C, f= const. gilt 

| für die gerade Kreiskegelschale mit in ihrer Symmetrieachse liegendem Momentenvektor. Die 
Bedingung des zweiten Ansatzes ce f= f” ist für f=sinAy-+ cosAy erfüllt®). Es muß nun die 
Kegelschale aufgesucht werden, bei der sich der Hebelarm s® der Schubspannung 7,, in bezug 
auf den Vektor des an der Kegelspitze eingeleiteten Momentes M mit cos Ay oder sin Ay ändert. 
Das ist wieder die gerade Kreiskegelschale und zwar mit zu seiner Symmetrieachse senkrechtem 
Momentenvektor. Es ist anscheinend nur in diesem Falle möglich, zwischen dem Hebelarm der 
Schubkraft, der beieiner mit ihrer Momentbelastung gegebenen Kegelschale eindeutig bestimmt Br: 
ist, und der Schubspannung unter Anwendung obiger Variation eine Beziehung abzuleiten. er 
Hierbei muß wegen Gl. (14) eine Achse des Koordinatensystems, in dessen Richtungen ein beliebig Er, 
gegebener Momentenvektor zu zerlegen ist, mit der Symmetrieachse des Kreiskegels zusammen- 2 
fallen. 

Während bei der Belastung der Kegelspitze durch Kräfte das angewandte Verfahren eine Sy 
eindeutige strenge Lösung des Spannungszustandes für jede beliebige Kegelschale angibt, führt a 
es bei Momentbelastung nur bei der geraden Kreiskegelschale zum Ziel. 

Zur Auffindung des Spannungszustandes bei der allgemeinen Kegelschale mit Moment- 

! belastung’ ihrer Spitze könnte als Näherungsverfahren für die Schubspannung ein den Rand- 
bedingungen genügender Ansatz gemacht werden, etwa in Form einer Fourierreihe, deren Kon- 
stanten so zu bestimmen sind, daß die Formänderungsarbeit zum Minimum wird. Dieses Ver- 
fahren ist nur zu empfehlen, wenn mit Sicherheit eine gleichmäßige Schubspannungsverteilung 
über die Schalenstärke zu erwarten ist. Wird z. B. eine auf Torsion beanspruchte Kegelschale 
längs einer Erzeugenden aufgeschnitten, so bildet sich eine Spannungsverteilung aus, bei der die 
Schubspannung nicht mehr über die Stärke des Kegelmantels konstant ist, sondern ganz im 
Gegenteil, sie wechselt sogar ihre Richtung. Dann ist obiges Berechnungsverfahren unbrauchbar 
und es müssen andere Wege beschritten werden, auf die hier nicht näher eingegangen werden soll. 


2) Die Hyperbelfunktion scheidet aus, da für sie keine stetig in sich zurücklaufende Kegelform existiert. 
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$4. Die Ei ptiaee Kegelschale bei Belastung ihrer Spitze durch Kräfte 
Legt man den Koordinatenanfangspunkt eines #-y-z-Koordinatensystems in die Kegel- 
spitze und sorgt ferner wegen Gl. (12) dafür, daß die «-z-Ebene und die y-z-Ebene Symmetrie- 


ebenen des Kegels sind, so lautet die Gleichung dieses Kegels (s, Bild 2): ne : 
22 1 rn TE 
—— Sr BES | 
An RR 


Sale) 
ßı und ß, sind die Winkel, die an der klei- 
nen und großen Halbachse der Ellipse von 
der Erzeugenden und der z-Achse einge- 
schlossen werden. Schreibt man Gl. (19) 
in Zylinderkoordinaten 2,r,p um, und be- 
zeichnet den von. einer Erzeugenden des 
Kegels und der z-Achse eingeschlossenen 
Winkel mit ß, so folgt (Bild 2): 


Bild 3. Die geometrische 22 tg? ß 
Beziehung zwischen den Winkeln tg? ß — 2 (20) 
dy,dy und dß 15 cos? .. ar 


ge! 2 
worin e® = nn die numerische Exzentrizität bedeutet. Es 
a 
ist oft leichter, den Spannungszustand als Funktion von g als vom 
Bild2. Elliptische Kegelschate Winkel y auszudrücken. Zwischen diesen Winkeln besteht, wie man 
: aus dem rechtwinkligen Dreieck ABC (Bild3) abliest, bei der all- 
gemeinen Kegelschale die Beziehung: 


(SUB Ir dr erde et Be (21). 


Eine einfache Umformung der Gl. (21) ergibt mit r = s- sin ß in Anwendung auf die ellipti 
Kegelschale nach Gl. (20) 


m we 
u 1—2co2p+tg?f, "(1-8 p+tgß)? (I —zog) f 


Betrachtet man zunächst eine Last P,in Richtung der 2-Achse, so wird der von einer Erzeugenden 
und der Last eingeschlossene Winkel £, identisch mit $. Der Spannungszustand os.) ist dann M 
nach Einsetzen von cosZ, als Funktion von durch Gl. (11) eindeutig bestimmt hi 


cos. 
=, Fr - ’ | 
fra zes 1 
9-0 x j 


Unter Anwendung von Gl. (20) wird cosZ, durch @ ausgedrückt: 


@ 
ER Verer 
8\z s6 28 Bew 
EN tg? ß, tg? B,etsin?pcos?p ; 
| Hier. Ve ee ee. 
Yen 


mit a=1—.e2cos?y. 


„. Für =, wird e=0, die elliptische Kegelschale geht in eine Kreiskegelschale über vom 
Öffnungswinkel 2, und Gl. (23) vereinfacht sich zu s = 


ung cos 1 
Ole) = — : SERBIEN, =— eng: SE er (24), 
f «os? Bsin Bd wenan 
go=( 


wie man auch auf einfacherem Wege finden-kann. 


1 Denkt ı man Sieh a Eee Last Pr, e die. Re Bet She liest aus sBild 2 die Beziehung. | 
co - 2: sin 2 cos IF ab, so folgt der. Spannungsverlauf aus. | 


ta f 


IR 


a 


Rp cos @ 


+2 h BR h | d a t 2 [92 B,. # N 4 u i% - y | ii Br t A 
B- ! - y Data) er sd re . hr . oo. a (25)... & ta 
AN HE EEE RR ig, Dr en Be: - 
AUS EN ve „erteg? & riet? 2“ | 0 

E . 2 Für die Kreickepeldchäle a sich GL. (25),zu‘. ve Da | % EN R I n” jr 
92,8 Er P NER RI en cosp ne ea 
Be RL 3 en) rt REN er sin? ß A ANE (26), j Rs 

Es ae n sin’ß ost Ei : / N 

F A en : Bet; ie 
2 - wie auch schon. bekannt ist, Es ist Rberfidssig, Ken ar für eine Last i in Richtung N: 


der anderen Halbachse der Ellipse auch noch anzuschreiben, da man ihn aus Gl. (25) durch 
'Vertauschen von ß, mit ß, leicht findet. fee TR IE 
Es darf dabei nicht übersehen werden, 


$% P x $, 
4 daß sich damit auch & ändert. Bild 4 5% 
E zeigt für eine elliptische Kegelschale A 
R vom Achsenverhältnis 2:1 und den , 
Winkeln tgß,=0,5; tgß,=1,0 den ws 
E Spannungsverlauf o, für die Last P,, te 
0 & .» El { X 
Bert rund P; in dimensionsloser ‚Dar- BB.’ 
a . r . er f 1 tar 
f stellung. Aus Symmetriegründen ge- Be 
3 nügt es, die Spannung nur über Y,der RE 
# Abwickelung der Schale vom Keilwin- BERN. BI 
kel 2, entsprechend den Winkeln 9=0 Er neuerer ihr oin 
x . bis 9% anzugeben, und zwar als Viertel dor. a 
2 Iptischen Kegelschalen 1 "Ri ER, 
Funktion von y, da die Darstellung o, om) Aöhsennrhiftnie 8-1 und A Fe 
| als Funktion von @ bei großer Exzen- en let Br einer  elliptischen. Kogelschale Ad \ 
N aa nach Bi ält- R “ 
. trizität e ein zu verzerrtes Bild gibt. P, und P, 2610 51 und den ee 3 
8Pı =0,1, tg, =1,0 unter . re 
Der Spannungsverlauf o, längs eines den Tasten ma PL Die R 
Schnittes s—= const. war in dieser Weise zu erwarten und Dt Spannungen für einen Keil vom EN 
keine Überraschungen, NE Ss 
Behält man den Öffnungswinkel ß, —- 45° an der großen Halb- flachen, Ellipse, wurden ge- Ra 
‘ striche . 
| achse bei, so daß tg ß,— 1,0 und verkleinert ß, zutgß,=0,1,so NE 
. erhält man eine flache Kegelschale vom Achsenverhältnis 10:1. Für eine Last P, parallel zur 


großen Halbachse und P, in Richtung der Mittellinie der Kegelschale folgt ein Spannungsver-, 
lauf o, nach Bild 5. Mit weiter zunehmender Verflachung der ER wird endlich = —=0 
und der Spannungsverlauf geht in denjenigen eines Keiles 2 5: 
[4, S.7] vom Keilwinkel 2,=% über, wie in Bild 5 ge- | 
strichelt angegeben wurde. Wie sich leicht nachweisen läßt, 
kann die auf anderem Wege bereits gefundene Spannungs- 
verteilung in der unendlichen Halbebene und im Keil [5, 
S. 294ff.] auch unter Anwendung von Gl. (11) als Grenzfall 
einer flachen elliptischen Kegelschale abgeleitet werden. 


3ER r 4 
ei 


‘ 5. Die offene Kegelschale bei Belastung ihrer 
Spitze durch Kräfte 
In diesem Beispiel soll die Anwendung der Gl. (11) und 

(12),auf eine offene Kegelschale gezeigt werden. Um das 
Wesentliche herauszuschälen wird, der einfache-Fall einer 
halben geraden Kreiskegelschale nach Bild 6 behandelt. 
Die Spannungsverteilung läßt sich für jeden beliebigen 
Kraftangriff in der Kegelspitze nach dem in Abschnitt I, 
$2 Gesagten sofort eindeutig unter Anwendung von Gl. (11) 
lösen, wenn die Lage des Achsendreikantes bekannt ist. Die e 
Richtung der x-Achse findet man in diesem einfachen Bei- nn TEE 
spiel sofort als senkrecht zur Symmetrieebene. Zur Bestim- gerade Kosfokagemenalk GERNE, 


xX;X 


bekannter Weise zu superponieren. Liegt der Momentenvektor in der Symmetrieachse des Kegels, a 
so kann man die Spannungen sofort anschreiben: ii 
| fi 
a TE MO | 
wTgahr 2d-2nsinip © Genla 7 DER WATER 028). H 
{ 


wo 


184 Sonntag, Beanspruchung der geschlossenen sowie auch offenen Kegelschale u 


mung der Richtung ‚der y- und z-Achse liest man unter Anwendung des 

Kugeldreiecke aus Bild 6 ab: | 

| cos £,= cos ßcos & +sin Psin & cos 
cos £, = cos ßsin & — sin P cos & cos 9, 


womit nach Gl. (12) mit dy - dp folgt: 


N Se 
| cos cos: dp= 0 = —sin2ß cos 24 +7008° Bsin25— 7 sin® sin 28 
e.- | | 
und damit dg22= dB — 7 1gß An re REN in ; 


Für = 45° folgt aus Gl. (27) &= 26°. 

Mit der eindeutigen Bestimmung der Lage des Achsenkreuzes hat man nun die Richtungen 
gefunden, in die die an der Kegelspitze angreifende Last zu’ zerlegen ist, um den Spannungs- 
verlauf nach Gl. (11) berechnen zu können, worauf hier verzichtet werden soll, um Wieder- 
holungen zu vermeiden. 


86. Die gerade Kreiskegelschale mit Momentbelastung ihrer Spitze 
Bei der geraden Kreiskegelschale vom Kegelwinkel 2 $ = const. sind aus Symmetrie- 
gründen nur zwei Belastungsfälle zu unterscheiden. Im ersten Fall liegt der Momentenvektor 
in der Symmetrieachse, im zweiten Fall steht er senkrecht zu ihr und geht nach Voraussetzung 
durch die Kegelspitze. Greift das Moment für einen gegebenen Belastungsfall in anderer Weise 
ansder Spitze an, so ist es in obige zwei Richtungen zu zerlegen und die Spannungen sind in 


rn Reg 


Steht der Momentenvektör senkrecht auf der Symmetrieachse des Kreiskegels, so ist der 
Spannungszustand nach $3 Gl. (18) durch den Ansatz f, = ec, cosAy = ® eindeutig bestimmt. 
Betrachtet man Bild 2, worin entsprechend einem geraden Kreiskegel $,—ß, zu setzen ist, und | 
liegt der Momentenvektor in der &-Achse, so ist das Moment der Schubspannung in bezug auf 5 

| 


diese Achse 
dM=rsy:Öcospscosß-sdy mit cosßcosp=6. 


Da ß = const. ist, gilt = En y n y und es folgen mit = ©, = cos ß cos ud y die 
Spannungen f 4 gr 


IT Re 
mM 2.008 — y Mm 2sin —Yy 
won sin2ß” ZT ed sm2PBsinß ey, Wer r 29). 


$7. Die elliptische Kegelschale mit Momentbelastung ihrer Spitze 
Es soll in diesem Beispiel das in $3 angegebene Näherungsverfahren angewandt werden. 


Eine elliptische Kegelschale nach Bild 2 wird durch ein Torsionsmoment belastet und zwar liegt 
der Moumentenvektor in der 2-Achse. Der einfache Ansatz 


1 


OT, SIE NIE AR 

Er; Tyan 2 sin —2p 
mit y=0 an der kleinen Achse der Ellipse genügt allen Randbedingungen und er wird für eine 
elliptische Kegelschale mit nicht zu großer Excentrizität eine befriedigende Näherungslösung 
liefern, anderenfalls ist die Hinzunahme weiterer Glieder der Reihe erforderlich. Die hier durch- 
geführte Rechnung bezieht sich auf den einfachen Fall einer elliptischen Kegelschale die nicht 
zu sehr von einer geraden Kreiskegelschale abweicht, etwa bis zu einem Achsenverhältnis 2:1. 


en Zul Da Ve De a au un 


a a EN ERS DU na sn, m a 2 25 
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ee BE ER 
‚in üblicher Weise aus der ;iimalbeöingung der Formänderungrarbeit. zu besti 
zur Vereinfachung die Poissonsche Konstante - = 0 gesetzt, womit also die 


vr 


{ 


| en vernachlässigt wird, so ist = 9 G, ind die Ointache Ditferentiation er Forman 
% ‚arbeit nach der Konstanten % ergibt nach einigen «Umformungen I 


cos = 2yp- a ; 


an 


Bei der ee Kreiskegelschale ist Ö, Rt damit wie das Integral, n = & Day .® EN a 


null, und c, verschwindet, wie zu erwarten war. Das Moment der Schubspannung Toy im bezug } Re 


auf die 2-Achse liest man aus Bild 3 ab: Baer ca 
dM = Ty öcosysin ßs?dy, "wobei. or D, sin? 07 BR 
| dy’ 


Für eine elliptische Kegelschale vom Achsenverhältnis 2:1 und den Kegelwinkeln ” Bis ==. * % 
tg ß, = 1,0 ergibt die Ausrechnung der Gl. (32) «= 0,12, und es zeigt sich nach Gl. (31), d ö 


die maximale Schubspannung ,, an der kleinen Halbachse der Ellipse auftritt, und zwar i ne % 
sie bei unserem Beispiel um 27% größer als an der großen ‚Halbachse. Die größe Beanspruchung 


nach Maßgabe der Hauptschubspannung 7z = Y 05/4 ii, v findet man an einem der Punkte, für 
die dry/dy verschwindet _ 
Ri dTg _ b, 4Acy [A?c, sin A2%Y- cos 22y—(1 +c, cos A2Yy) sin A2y] 
dp 208 Valat+ z, 
Es folgt 


ET aA ee 


Welche Werte dee SAL HDE der Gl. (35) ein Maxinjurn liefern, findet man Heike Se Ein- 


setzen in Tg. ' 
Eine maximale Ei ortelihe liegt für c, [(& 2) — | < lan der kleinen Achse, y=Ounda 


2 
vor;ist der Wert > 1, so findet man den Ort von Tymax aus cos = v= 1/6, (2) rd | £ 
.1\& 


III. Zusammenfassung ' 


Bei der allgemeinen Kegelschale mit Belastung ihrer Spitze durch Kräfte treten nur Längs- 
spannungen auf, die in verhältnismäßig einfacher Weise eindeutig bestimmt werden können. Das 
Rechnungsverfahren kann ohne Einschränkung auf eine offene Kegelschale, d.h. bei der keine 
in sich zurücklaufende Schnittkurve s= const. existiert, übertragen werden. Besondere Be- 
deutung kommt dieser Erkenntnis in Anwendung auf eine Kegelschale mit eingeschnittener 
Öffnung zu, da die die Öffnung überstreichende Mantelstrahlen im em näherungsweise 
als spannungsfrei angesehen werden können. 

Bei Momentbelastung der Kegelspitze sind die Verhältnisse viel komplizierter. Eine ein- 
deutige Lösung in geschlossener Foım war nur bei der geraden Kreiskegelschale möglich. Für 
anders gestaltete Kegelschalen wird ein in vielen Fällen befriedigendes Näherungsverfahren an- 
gegeben. Verschiedene Beispiele zeigen die Anwendung der Lösungsverfahren. 
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Der Schlupffehler des Reibradintegrators. 

Bei einem für das Max-Planck-Institut für Physik 
gebauten Reibradintegrator untersuchte der Verfasser 
eine Erscheinung, die zwar auf einem einfachen 
physikalischen Zusammenhang beruht, aber wenig 
bekannt ist. Die Rolle der in Bild 1 gezeigten An- 
ordnung bestand aus zähhartem Stahl, die Scheibe 


Bild 1 


aus einem gewöhnlichen Maschinenbaustahl. Der 
Rand der Rolle war mit einem Radius von 0,5 mm 
abgerundet und poliert, %die Scheibe geschliffen 
und poliert; beide waren leicht gefettet. Die zy- 
lindrischen Lagerzapfen der Rolle waren poliert und 
liefen in Bronzebuchsen. Von der Rolle wurde ein 
wenig Kraft verbrauchendes Zählwerk mit schleichen- 
den Ziffernscheiben angetrieben. Der Anpreßdruck 
der Rolle war ca. 400 000 dyn. 

Beim Betrieb des fraglichen Gerätes zeigten sich 
unerwartet hohe Fehler, z.B. beim Planimetrieren 
eines rechtwinkligen Dreiecks mit achsparallelen 
Katheten. Wenn in einem solchen Fall die Rolle zur 
Verfolgung der Hypothenuse radial um 7,89 cm = 
Y,— 4, unter gleichzeitiger Drehung der Scheibe um 
10,46 rad = 2,— x, verschoben wird, während die 
beiden Katheten durch Verschieben des Rollensupports 
bei ruhender Scheibe bzw. durch Drehen der Scheibe 
bei ruhendem Rollensupport gebildet werden, so 
beobachtet man an Stelle der Rollendrehung 


Hua —an) 9 9104 
PA 


einen um etwa 0,38 rad — beinahe 3% ! — zu kleinen 
Wert. Der Fehler ist gut reproduzierbar und ziemlich 
unabhängig davon, an welcher Stelle der Scheibe der 
Versuch ausgeführt wird. Er verschwindet völlig, 
wenn gleichzeitiges Verschieben des Rollensupports und 
Drehen der Reibscheibe vermieden werden, " die 


endlichen Momentes M muß diese Rejbungskraft eine 

iale Komponente M/r haben, wo r der Rollen- 
radius ist. Daher weicht die obige Reibungskraft um 
einen kleinen Winkel 


in der. gezeichneten Richt I 5 
rüh punktes ab. In dieser Richtung tritt zwischen 
Rolle und Scheibe reines Gleiten ein, und 


AC 1 ae 
= I — —ydy? 2 (C0E &— 
dz = ah dy? + (y de)? (coe a—sinatanß) 
erteilt wird. Mit cosa = ydz/Ydy®+(ydz)* und 


sina = dy/Ydy? + (ydx)? folgt 


de= 1 (ydz —dytanß). 


Das erste Glied in der Klammer liefert das zu messende 
Integral, das zweite gibt nach Integration innerhalb 
der Grenzen y, und y, den Fehler 


_ 


M 
0) = NY) ag 


Löst man diese Beziehung mit den im Zahlenbeispiel 
schon gemachten Angaben sowie mit der Annahme 
e=0,1l nach M auf, so erhält man den für das be- 
trachtete Gerät sehr plausibeln Wert M=11 020 dynem. 


Zu beachten ist, daß der Drehwinkel x der Scheibeim . 


Ausdruck für den absoluten Fehler nicht vorkommt. 
Der relative Fehler der Integration kann also durch 
Vergrößerung von x gesenkt werden. Durch Plani- 
metrieren eines aus der Diagonalen und zwei Seiten 
der Schreibfläche bestehenden Dreiecks kann der 
Reibradintegrator auf seinen Schlupffehler geprüft 
werden. 


Göttingen. K.Beyerle. 


Zur Praxis der numerischen Quadratur. 

Die Korrekturglieder für die Mittelwertformeln der 
numerischen Quadratur eignen sich meist wenig dafür, 
beim praktischen Rechnen den Fehler des Ergebnisses 
zu bestimmen, Ihren deutlichen Ausdruck findet diese 
Tatsache in den Beispielen ‘der Lehrbücher, deren 
Genauigkeit durch Vergleich mit dem bekannten Wert 
festgestellt wird. Diese Schwierigkeit läßt sich nun 
dadurch umgehen, daß man Einschließungsformeln 
benutzt, Aus der großen Fülle der bekannten Qua- 
draturformeln findet man die hierfür geeigneten leicht 
durch Beachtung der Vorzeichen und Koeffizienten 
der Korrekturglieder. Für die Fehlerbestimmung bei 
Anwendung der Quadraturformeln mit Differenzen 
hat Charlier [1] eine Regel aufgestellt. 
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ETEREET 


EHYt Aus w 


erhalten wir 


Ne ‚terra 
Tante die Simpson sche Formel 


ran ent var 


es Ba CIE 


Als Bigsnsugpforme eignet, sich Be a x 
2 ie x 1-77 Sr “Blut sudt und 


N 
rn Br 6 a) 10) (n). 
‚Als Br Boll aa 


Sie are 
ee 


berechnet werden. A ‚den Ordinaten EN 


Yo il ‚0000 
Yıra = 0,6667 


Yıc = 0,5454 


0,6945 ) In 2 ) 0,6926 


und da ihre Fehlerbeträge sich annähernd wie 


3:1 
verhalten, ist: ji drei Dezimalstellen genau n2—= 


0,693. 


H. Nach Durasıal Peano [3] ist 


Tao wor +5 YA 


eo 
300 7! 


®. 
od, 5 


wo BEER ER 
‘a = 0,27639 32023, 


B = 0,72360 67977. 
Als. ah pn eignet sich die Gauss es 


ie I: +8 Ya 


% 2 
em, 
wo 
y = 0,11270 16654, 
ö = 0,88729 83346. 
Als Beispiel sei wieder In2 berechnet. Mit den 
Ordinaten 
Yo = 1.000000 
Yı = 0,500000 _Yıra = 0,666667 
Yu 0,183458 y, = 0,898713 
Yg= 0,580179 y; = 0,529858 


erhalten wir 


0,693182 ) In2 ) 0,693122 


und da sich ihre Fehler verhalten wie 4:3, ist auf 
fünf Dezimalstellen genau In 2 = 0,69315. 


III. Die nächsthöhere Turazza-Peano- For- 
mel are 


et 


+yıl+ 64 yıla + 49 [y, + Yo 


1 (b—a)? 
3920 : 91 


70 
HOIGE 


} sche: 


- sparen. 


‚ einen Fehler, der von 6, Ordnung klein ist. Hieraus RR 


Er 82 0,17207 31047, 
= 0,82732 68353. 


Als RER eignet sich wieder e eine 


I 0a) {4 Int 9 ” Bi, Br 223 i ; = re 
2900 RE 2 . 
wo u s4 


A = 0,17392 74226, B= 0 ‚32607 25T, 
“x = 0,06943 18442, m .0,33000 94782, 
A EU) ‚93056 81558, v= 0 ‚66999 05218. % 
Als u soll ee Er 
: RE Ed 
de 
Az 
EN I x 
berechnet Werden: Mit ae Ordinaten. el) 


”. 
6 


“ = 1,00000 000 a 314 25 
Yı =1,15470094 = 112973745 
un 03279556 er RE, 
„ = 1,00374 797 Y „= hasse. ee 
n = 1,09838 222 —=10612237 00000 
erhalten wir ° = Be © 


'0,52359 889 x rn ax 0, 52359 869 


und da sich ihre Fehler Et wie 5:4 verhalten, ist e 
auf 7 Stellen Fe 


= = 0,52359 88. 

Diese Einschließungsformeln werden für die zieikten 
praktischen Bedürfnisse genügen. Höhere Genauigkeit. 
läßt sich stets durch Intervallteilung erreichen. Den- 
noch verdient die Frage gestellt zu werden, ob es nicht 
möglich ist, den Formelfehler ohne Berechnung neuer 
Ordinatenzufinden und die durch die irrationalen | 
Abszissen oft notwendigen Tafelinterpolationen zu 
Solche Formeln sollen jetzt für 5 und 11 
Ordinaten angegeben werden. 

IV. Bilden wir vier äquidistante Teilintervalle, so 
ergibt die doppelte Anwendung der Sim psonschen 
Formel einen Fehler, der von 4. Ordnung klein ist. 
Dagegen hat die Co t essche Formel für 5 Ordinaten 


folgt, daß wir mit ausreichender Genauigkeit den 
Fehler der Simpsonschen Formel durch ihre 8 
Differenz mit der Cotesschen bestimmen können. TIER 
So entsteht mit der Abkürzung E FE 


w=1la4 129) | 

das folgende Formelpaar _ ; 
2 a . 

I= a 


ae SSETIN td —ıln+yml+6yd:. 


Das Korrekturglied K dient zur Abrundung der letzten 
Stelle [4]. 


Als Beispiel soll 


al 
Fr& 2) 


2 12 
0 
berechnet werden. Mit den Ordinaten 
Yo = 1,00000 Y; = 0,74616 


Yyı = 0,89260 
Y, = 0,81094 


Yı= 0,69315 


- 


ln well 


er A Pan:! S A, 
en N Br £ A 
» ’ x , 
188 Ent i ö RN s m 


vier Stellen genau 


ne 
V. Die numerische Quadratur mit 10 


une Ran 
o = b— a/l0 ist besonders TER da sich am 
und Interpolationen spart. Mit der Abkürzung 


vi=f(a+io) 


lautet die Cotessche Formel für 5 Teilintervalle [5 


a+5w 


EIL: 


= (19 yo+ 75 y1 + 509.4 509, oyı+ 104) 


I a {yD — yD) 


12096 
und hiermit folgt 
TE» 


576 


R 55 
} 


19Wy+ Yı)+ dlyı + y+ Yet Yyol 
+50[y + Yy3+ yr + Ysl + 38 4} 


EI b—a 6 y® 
12096 \ 10 w 


Anderseits ist nach Sheppard [6] 


2 2 Das 


NN 
"1512 


Setzen wir in leichtverständlicher Abkürzung 


IT 8lyo + Yıo) + 35 Wyı + Ya + Yr + 95] 
+l1ö[y +y4+Ys+ Ysl + 36 9} 


[b—a\°y, 0 
( 10 ) 70.08 


I=C—EBc=8S—Es, 


so ist 


Bo Es= (08 


und wegen 


Eo :Es =45 :33 


finden wir 


33 
Es = 75(0—8). 


Diese Fehlerformel dient wieder zur Abrundung der 


letzten Stelle. 


Da von den Korrekturformeln nur das erste Glied 
berücksichtigt wurde, ist die Gültigkeit dieser und 
der unter IV angegebenen Fehlerbestimmung davon 
abhängig, daß b—a hinreichend klein ist; diese 
Forderung kann aber durch Teilung des Integra- 
tionsintervalls stets leicht erreicht: werden. 


Als Beispiel soll 
200 000 


100.000 


berechnet werden. (Auf dieses Integral hat bereits 
Gaussseine Methode angewandt.) Mit den Ordinaten 


Ye = 8345,20502 


Yo = 8685,88964 
y, = 8614,57357 
Y. = 8550,48207 
Ys = 8492,36004 
ya = 8439,24751 
Ys = 8390,39461 


erhalten wir 
C = 84.06,24321 


2 
=, du 


5inl0 + Inu 


Yz = 8303,19708 
Ys == 8263,97644 
Yo = 8227,21640 
Yıo = 8192,64336 


5 = 84.06,24319 


Es = 5.107, 


Daher ist bis zur letzten Stelle genau / = 8406,2431. 


82251 —0,00002, d.h. es ist auf 


| len! 
Turazza, . Ist. Veneto 5 (1855), p. 255 bis 
En Peonor Attı di Tora 27 II E 
p. 608—612. (Peano hat die zugehörigen 2 
nome abgeleitet.) 
[4] Mit dem Rechenaufwand, den die vonv.Sanden 
rear Analysis) empfohlene rohe 
immung der Simpsonschen Formel ver- 
ursacht, läßt sich auf diese Weise eine weit bessere 


[5] Dies Hosted wurde von A. Klingst, 
Ber. Akad. Wien, Abt.I, 150 (1941), 8. 221 bis 


256. 

[6] Diese Formel findet sich in anderer Schreibweise 
und ohne Restglied zuerst bei Sheppard 
Proe. Lond. Math, Soc. 32 (1900), p. 258—277. — 
Shovelton, nach dem Whittaker und 
Robinson (Caleulus of Observations) diese 
Formel benennen, hat sie erst 1914 in der gleichen 
Zeitschrift (ohne Restglied) abgeleitet. 
Göttingen. F.K.Rubbert. 


Der mittlere Fehler bei der Mittelwertbildung 
aus Wertgruppen. 


o 
In der Fehlerausgleichsrechn! lt für die Berech- 


nung des Mittelwertes x von bnissen x; mit 
den jeweiligen Gewichten p; die Gleichung x 
| [p: x] 
4 i 
ST ERR m 


und für die Berechnung seines m. F. (mittleren Fehlers) 
M die Gleichung 


M= [p: v5] (2) 
} (a—1) [pi] ’ ’ 


worin [....] die in der Ausgleichsrechnung übliche 
Schreibweise für das Summenzeichen £.,... darstellt 
und v; = #— x; die Verbesserung des Meßergebnisses 
x; bedeutet. Diese Gleichungen und ebenso die folgen- 
den Betrachtungen gelten auch in der Statistik für das 
arithmetische ittel einer Reihe von Kollektivgegen- 
ständen und seine Streuung. 

Bei der Anwendung der obigen Gleichungen liegt 
häufig der Fall vor, daß die Ausgangswerte z; der Rech- 
nung nicht mehr Ergebnisse einzelner Messungen, son- 
dern selbst schon Mittelwerte von in Gruppen zusam- 
mengefaßten Einzelmessungen sind. Solche Gruppie- 
rungen von Meßwerten können aus verschiedenen Grün- 
den vorkommen, Beispielsweise kann der Fall vorlie- 
gen, daß Meßwerte des gleichen Gegenstandes von ver- 
schiedenen Stellen unabhängig voneinander ermittelt 
wurden, die man nun zusammenfassen will. In solchen 
Fällen kennt man aus den-Mitteilungen der verschie- 
denen Beobachter die Einzelwerte häufig gar nicht 
mehr, sondern nur noch ihre Mittelwerte und deren m. 
F., vielleicht noch ihre Anzahl. Oder es soll ein sehr 
umfangreiches Material durch mehrere Bearbeiter aus- 
gearbeitet werden; man überläßt dann jedem von ihnen 
eine Gruppe und benutzt weiterhin nur noch die Grup- 
penmittelwerte und deren m. F,. Oder auch man redu- 


| eren 
sen verschiedenen mögli« n 
s Gruppen- festen Wert M als Mittel 

nit e d seinen m. können erheblie ‚ausf 
erhalten, als wenn man Fer den 
gegangen wäre. 
"bekannt, daß man für den Mittel. 
en) dasselbe Ergebnis bekommt, ob Berechnung zu e 


Gl. (1) unmittelbar aus. den Einzelmes- Es ist nun möglich, diesen ER ER Wert M 
et, oder ob man zunächst Gruppen- des m.F. mit Hilfe einer neuen Gleichung ı auch ohne 


ch Gl. (1) und erst aus diesen unter Be-  Kenntnii 
s der Einzelmeßwerte aus G: Be 1 
ler Gewichte einen Gesamtmittelwert ten zu berechnen, wenn wenigstens AN ui 


gste, da seine re ı 
von Ausgangswerten er Igt. Er 


bildet. : | 
00 Bsist aber weniger et daß die dl. @) für den Be Eon De Be 


x m.F. des Gesamtmittelwertes bei der Berechnung aus ch 
Eu Gruppenmittelwerten ein abweichendes Ergebnis ge- y ISR ee 
...  genüber der en ng aus. den Einzelmeßwerten 
liefert. Würde man mit denselben Meßwerten eine neue —— 

SE u lue vornehmen, so würde man bei glei- RR (ul a [pl a 
E: , : he ) nie BER a heranie ; WA EB hist : 3% 
{ Ri} Berechnung aus Einzelmeßwerten et, Fr 
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Mur RL : Ast EEE 


[P: v7] in en, Pi ;M a 


N M= 


IE 
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: wert |: Gesamtmittelwert « sein m,F. M Gruppenmittelwerte &bisa, . ihre m Hr ‚M, bisM, ER ; 
£ e n.,Gl,(hjl, 5 = n. GM) n..Gl.(1) | Gl. (2) 


RE; z;!) D | pi | pizi vi | piv; | | | 
E.; N ee a a 111 | —44 | 194 RE ASt IP Sa Pre 
4 156 | 2 1 156 | —20,6 | 424,4 1 2 1 156 | +125 | 156,3 
3 | PR En ee 
B: 192 | 3 1 192.1 — 15,4 |. 237,2 | 1 192 7 1900.00. 807 
E 168.14. 9° 1 165 | +11,6 | 134,6 A 1 165 13.3. ,..1769 
Be - 178 | 5 1 ee 2,0 3 1 178 EEE EEN. 
E ae | 2. ern er 
3 157 | 6 1 157 +19,6 | 384,2 [4 1 157 +22,3 | 497,3 
3 178 | 7 1 1787-1. 54 RER 2 1 178 En ee 
2 180 | 8 14 180.210 3.4 11,6* ER Ne 1 180 07 0,5 
£ 202 | 9 1 202 | —25,4 | 645,2 74 1.4 209 —22,7 | 515,3 
: Zr Be, 1089 eo > | 1860,6 a a te: 
4 1589 337 
=, = 116,6 M | Er 1 re = >= 12,50 
r 8 F 
= 5,08 3 Se 
£ 535 364,7 
2 = 183 I0,=) 5; = 1,80 
| 717 = 
3 Ber =lds = 9,20 
; Berechnung aus Gruppenmittelwerten . 
| %i, Pi, nı aus III, M; aus IV S 
V vI S+ | vu ; 
Gruppen- 27 BE EEE TEE 
: Gesamtmittelwert x seinm. F, M sein m.F,.M 
a he nach GI. (1) nach GI. (2) nach GI. (3) : 
wi i |» | ps vi | pie | (ml) pi M} 5 
168,5 1 2 337 AL 131,2 2 12,50 312,6 
178,3 2 3 535 N 8,7 3 7,80 364,7 
179,3 3 4 TLF Ey 29,2 4 9,20 1014,8 
Being: 1889,87] el H1e9 17.5029 - 1692,1 
ar NM = nn. = 3,06 M= oa os 
9 29.9 
Ergebnis übereinstimmend mit I abweichend von II ° übereinstimmend mit mit Il 


2) Die Ausgangswerte stellen die Höhen von neunjährigen Kiefern in cm dar. Sie stammen aus H. Vater, Die Ausgleichsrechnung bei 
Bodenkulturversuchen, Mitt. a. d. sächsischen forstlichen Versuchsanstalt zu Tharandt, Bd.II (1925), S.8, "und wurden als Mittelwerte 
aus je 14 Messungen gebildet, die aber hierim Interesse einer gekürzten Darstellung wie Einzelmessungen betrachtet werden. 


TO 


' Sie\ist eine allgemeinere Gleichung, die die Gl. (2) 
als Sonderfall in ch schließt und in diese übergeht, 
wenn jede Gruppe nur aus einem Meßwert besteht, also 
die Gruppenstärken n; = 1 werden. 

Ihre Ableitung erfordert einen größeren Aufwand an 
Gleichungen und neuen Bezeichnungen, sie soll aus 
Platzgründen nur kurz angedeutet werden. Sie beruht 
darauf, daß bei der Berechnung des m. F. aus Einzel- 
meßwerten nach Gl. (2) nicht die Kenntnis jeder einzel- 
nen Verbesserung v;, sondern nur die Kenntnis der 
Summe [’p: v3]erforderlich ist. Diese Summe kann man 
auch aus den Verbesserungen der Gruppenmittelwerte 
und aus deren m. F. M; berechnen. Die dafür notwen- 
dige Gleichung gewinnt man, wenn man davon ausgeht, 
daß sich die Abweichung jedes einzelnen Meßwertes 
vom Gesamtmittelwert aus seiner Abweichung von 
seinem Gruppenmittelwert und dessen Abweichung 
vom Gesamtmittelwert zusammensetzt. Diese für je- 
den Meßwert gültige Gleichung führt durch Quadrieren, 
durch Multiplizieren mit dem jeweiligen Gewicht und 
durch Addieren über sämtliche Meßwerte zu der ge- 
suchten Beziehung zwischen den gewogenen Fehler- 
quadratsummen, die im Endergebnis die Gl. (3) liefert. 
Die Übereinstimmung mit dem aus Einzelmeßwerten 
berechneten m.F. ist exakt, dabei der Ableitung der 
G1.(3) keine Vernachlässigungen vorgenommen wur- 
den. Etwaige kleine Abweichungen beim Zahlenrech- 
nen beruhen nur auf Abrundungsfehlern. 

Die neue Gl. (3) gibt uns die Möglichkeit, den zuver- 


lässigsten Wert M für den m. E, des 
wertes zu berechnen 


Berechnung ai Es ist nicht zu empfehlen geh 
ee Gen leer Na 
Gl.(2) zu koghge (s. Zahlenbeispiel nnlenbersniel), Wert ur } 


Aufgabe keine siche- 
ven Anhaltspunkte für die Wahl der Gewichte p; der R 


Gruppenmittelwerte und für 


die yes Et 
besitzt, kann man in den Gl. (1) bis (3) Dan mer 


in üblicher Weise umgekehrt i 
ten der m.F. M; annehmen. In diesem Falle nimmt 
die G1.(3) die folgende, für die Zahlenrechnung beson- 
ders bequeme Gestalt an: 


I Ar 
» .(ia]—1) la 


worin n die Anzahl der Gru ittelwerte bedeutet. 
Hier ist noch die Gesamtzahl [n;] aller MeBwerte ent- 


M= (3a), 


halten, die man aber — falls nicht bekannt — nur. 


ößenordnungsmäßig richtig zu schätzen braucht, da 
ihr Zahlenwert das Ergebnis wenig beeinflußt. Nur 
wenn auch die m. F. M; der ig Fra un- 
bekannt sind, ist man weiterhin auf die Benutzung der 
Gl.(2) angewiesen. 


Speele (Kr. Hann. Münden. W.Weydanz. . 


BUCHBESPRECHUNGEN 


Prof. Dr. P. Riebesell, Mathematik des täg- 
lichen Lebens. 828. Hamburg 1948. Verlag 
Hammerich u. Lesser. 

Diese Plaudereien über Mathematik im täglichen 
Leben sind in den Jahren 1942—44 in 4 Auflagen er- 
schienen und in dieser Zeitschrift Bd.22 S.301 be- 
sprochen. In der neuen Ausgabe sind einige Abschnitte 
(Relativitätstheorie, Atomenergie) hinzugefügt. Auch 
diese Auflage wird ihre Leser finden. Für einen Neu- 
druck würde sich eine sorgfältige Durchsicht empfehlen, 
um den Ausdruck an manchen Stellen mehr zu präzi- 
sieren und Versehen zu beseitigen; z. B. wird von der 
Keilschrift der Agypter gesprochen, es wird gesagt, daß 
„im allgemeinen bei den Kartenprojektionen die Groß- 
kreise als gerade Linien erscheinen‘‘ usw. 


Dresden. Willers. 


Dr. Wilhelm Sehlink (Prof. a..d. Techn. Hochschule 
Darmstadt) und Dr.-Ing. habil. Heinrich Dietz (Doz. a. 
d. Techn. Hochschule Darmstadt), Technische 
Statik. Ein Lehrbuch zur Einführung ins technische, 
Denken. Vierte und fünfte Auflage. 431 8. m. 511 Abb. 
Berlin, Göttingen, Heidelberg 1948. Springer Verlag. 
Preis brosch. 27,60 DM. 

Das in neuer Doppelauflage vorliegende Werk ist seit 
seinem ersten Erscheinen in den Kreisen der Ingenieure 
und Studierenden so bekannt geworden, daß es einer 
neuerlichen Empfehlung entraten kann. Es enthält die 
grundlegenden mechanischen Fragestellungen der Bau- 
statik, soweit sie der Mechanik der starren Körper an- 
gehören, in einer didaktisch erprobten und vor allem 
den Studierenden der ersten Semester angepaßten Dar- 
stellung. Es werden mit zahlreichen Beispielen ebene 
und räumliche Kräftesysteme in Verbindung mit 
Gleichgewichtsaufgaben bei geraden und gekrümmten 
Stäben, Fachwerken und Rahmen, sowie bei Gemischt- 
systemen behandelt, wobei auch auf die Ermittlung der 
Einflußlinien für Auflagerkräfte und innere Beanspru- 
chungsgrößen. bei bewegten Lasten eingegangen wird. 
Bei symmetrischen Bauweisen wird auf die mit den 
Symmetrieeigenschaften verbundenen Vereinfachungen 
hingewiesen. Bei der Behandlung räumlicher Probleme 
wird auch die „zugeordnete ebene Abbildung‘‘ und ‚‚die 


duale Kräfteabbildung‘‘ bei beliebig zerstreuten Kräf- 
ten erklärt. Es mag jedoch auffallen, daß dabei die für 
die Raumstatik grundlegenden Begriffe, wie Dyname, 
Zentralachse u.a. nicht erwähnt werden. 

Mit. besonderer Betonung muß jedoch die Tendenz 
des Werkes hervorgehoben werden, zugleich mit dem 
sachlichen Inhalt der Statik auch eine Einführung in 
das Wesen der Konstruktion zu geben und so eine 
Grundlegung für die Ausbildung des technischen Den- 
kens dös werdenden Ingenieurs zu schaffen. Diese Ten- 
denz kommt schon im Titel des bestbekannten Werkes 
zum Ausdruck. 


Karlsruhe. Th. Pöschl. 


Dr.-Ing. Hayo Föppl, Messung der Span- 
nungen und der DIRSEInSERE Verfor- 
mungen an oberflächengedrückten 
Probekörpern, (Mitteilungen des Wöhler-Insti- 
tuts Braunschweig Heft 41.) 67 S. m. 40 Abb. Braun- 
ae 1948. Verlag Friedrich Vieweg u. Sohn, Preis 
4,— x 

Bekanntlich wird durch das Oberflächendrücken, das 
auf verschiedene Arten (durch Hämmern, Rollen, Ab- 
blasen mit Stahlgrieß u. dgl.) ausgeführt werden kann, 
die Dauerhaftigkeit eines Werkstücks wesentlich er- 
höht. Zur Ermittlung der in der Oberflächenschicht 
beim Drücken entstehenden Spannungen sind verschie- 
dene Wege angegeben worden (z. B. von ©. Föppl, 
G.Sachs, H. Chitz usw.). Die hier zur Überprü- 
fung dieser Methoden angestellten Versuche wurden in 
der Weise ausgeführt, daß die Probestücke in ihrem 
zylindrischen Teil geätzt und das entstandene plastische 
Verformungsgebiet stufenweise wieder abgenommen 
wurde, Durch Messung der Länge des Probekörpers 
und der Dicke der abgeätzten Schicht werden in jeder 
Stufe die für die Spannungsberechnung notwendigen 
Beobachtungswerte gewonnen. Durch Veränderung der 
Art des Drückens (Anpreßkraft, Vorschub, Abrun- 
dungsradius der Rollen und Zahl der Überrollungen) 
werden die jeweils günstigsten Ergebnisse für das Ober- 
flächendrücken ermittelt. 


Karlsruhe, Th. Pöschl, 


a ee 


Fi 


« > 
Br, 


ahn. Hoch. 


‚ Grundbegriffe, die zur Kennzeichnung des Spannungs- 
und Vera derertordemmieh 


lasteten Kreiszylinderschale mit allgemeiner Lösung für 
konstante, Wandstärke unter Flüssigkeitsdruck, die 
Membrantheorie der Zylinderschalen mit Anwendung 
auf die Zylinderschale mit waagerechter Achse unter 
Eigengewicht und auf Zylinderschalen allgemeiner 
Form, und schließlich Einzelheiten des Spannungszu- 
standes behandelt. Die zahlreichen Abbildungen tragen 
wesentlich zum Verständnis der formelmäßigen Ablei- 
tungen bei. | ! ; 
Vom Leser wird an Vorkenntnissen ungefähr das vor- 
ausgesetzt, was der Vorprüfung im Bauingenieurwesen 


LE 1948. 
sgeseben 
Wolfen. | 


an einer Technischen Hochschule entspricht. Da das 


Heft klar und leicht verständlich geschrieben ist, kann 
die Zielsetzung, dem Leser eine Einführung in das um- 
fangreiche Gebiet der Schalenstatik zu vermitteln, als 
gelungen bezeichnet werden. 


‘ Karlsruhe. a Th. Pöschl. 


Dr. phil. H, von Sanden (o. Prof. an der Technischen 
Hochschule Hannover), Praktische Mathematik. 
Teubners Mathematische Leitfäden, Band 44. 100 8. 
mit 17 Abb. Leipzig 1948, B. 6. Teubner, Verlags- 
gesellschaft. Preis 3,30 DM. 


Das Buch ist aus einer Vorlesung über praktische 
Mathematik hervorgegangen, die an der Technischen 
Hochschule Hannover neben den Anfängervorlesungen 
über höhere Mathematik -gehalten wird (eine sehr 
beachtenswerte Regelung!). Es bringt für den. Stu- 
denten der Ingenieurwissenschaften und den in der 
Praxis stehenden Ingenieur in fünf Abschnitten 
(Zeichnerische Behandlung von Funktionen; Der Satz 
von Taylor, Näherungsformeln; Integration, Differen- 
tiation und Interpolation; Statistik; Ausgleichsrech- 
nung) einen Querschnitt: durch die gebräuchlichsten 
‚Methoden der praktischen Mathematik. Das Büchlein 


setzt nur geringe Vorkenntnisse in höherer Mathematik 


voraus und ist klar und leicht verständlich geschrieben. 
Eine große Zahl von Beispielen stellen die nötige 
Verbindung zwischen den Fragestellungen der Praxis 
und den gebotenen Lösungsmethoden her. Besonderer 
, Wert wird dabei auf planvolle Anlage jeder Rechnung 
gelegt, und der ganze Rechnungsgang wird auf die 
gewünschte oder erzielbare Genauigkeit des End- 
resultates . ausgerichtet. Zu kurz erscheinen dem 
Referenten die Auflösung von Gleichungen (mit 2 Seiten) 
und der Anhang über harmonische Analyse behandelt, 
Man kann das Buch besonders auch den Studenten 
‘der Technischen Hochschulen empfehlen, an denen 
noch keine Vorlesungen über praktische Mathematik 


, für die technischen Fachrichtungen abgehalten werden. 


Dresden. Opitz. 
Dr. Wolfgang Haack (vorm. o. Prof. a. d. Techn. 
Hochschule Berlin-Charlottenburg), Differential- 
Geometrie. Teil II (Bücher der Mathematik und 
Naturwissenschaften. Herausgeg. von Dr. Henry 


v 3 R ” . 
olfenbüttel-Hannover‘ 


_ Wolfenb 
brosch. 7,—DM. 


Lie, Riemann, Levi-Civita usf. gekennzeichnet. 
Aus der Fülle ihrer Gedanken und Untersuchungen 
bietet der Verfasser in dem vorliegenden zweiten Band 
seiner Differential-Geometrie eine wohlerwogene Aus- 
wahl dar. Der Inhalt des nur 130 Seiten starken Bänd- 


chens umfaßt eine Einführung in den Ricei-Kalkül, _ 
eine Darstellung der invarianten Ableitungen der 
. Flächentheorie und in vier weiteren Kapiteln einen 
Abriß der differentiellen Liniengeometrie (Elemente, 


Regelflächen, Strahlensysteme und -komplexe). Eine 


‚ fruchtbare Lektüre dieses zweiten Bandes setzt die 
'Vertrautheit mit dem ersten voraus; aber auch. wenn 


diese Bedingung erfüllt ist, wird der Leser bisweilen 
nichts zu lachen haben. Das liegt aber dann weder. 
an ihm selbst — wie ihm zum Trost gesagt sei—, noch 


am Verfasser sondern allein an den unvermeidliichen 


Schwierigkeiten, die Fülle des Stoffes mit der Enge 
des Raumes in Einklang zu bringen. Das Mögliche 
in dieser Hinsicht hat der Verfasser getan und damit 
sein pädagogisches Geschick und seine überlegene 


Sachkenntnis erneut bewiesen, Die äußere Ausstattung 


des Notdruckes ist zeitgemäß bescheiden. 
. Freiberg (Sa.). Grüß. 

Prof. Dr. Erich Schönhardt, Vektor-Rech- 
nung-mit jeeinem Anhang über Ten- 
soren und über komplexe Zahlen und 
Zeiger. Herausgegeben von Dipl.-Ing. Werther 
Pavel(2. Aufl.). 295 8. m. 150 Abb. Stuttgart 1948, 
Preis brosch. 12,— DM. i 


Die. vorliegende Einführung in die Vektorrechnung 
ist in der Hauptsache für Ingenieure und Physiker ge- 
schrieben, daher wurde an einigen Stellen nur der Ge- 
dankengang der Beweise ohne eingehendere Ausführun- 
gen gegeben. Das Buch zeichnet sich durch eine große 
Zahl (155) guter Beispiele aus, zu denen meist ausführ- 
liche Lösungen gegeben werden. Ein Student, der 
möglichst selbständig diese Beispiele durcharbeitet, 
dürfte dadurch mit der Handhabung der Vektorrech- 
nung vertraut werden. Das Buch gliedert sich in einen 
verhältnismäßig umfangreichen Abschnitt über Vektor- 
algebra (etwa 80 Seiten) und einen über Vektorana- 
lysis (125 Seiten). Ein erster Anhang (30 Seiten) bringt 
die Grundbegriffe der Tensorrechnung, ein zweiter 
(40 Seiten) das Wichtigste über komplexe Zahlen und 
Zeiger. Ein alphabetisches Register erleichtert das 
Nachschlagen. 

Das treffliche Buch, das sehr sauber geschrieben und 
durch ein Flachdruckverfahren vervielfältigt ist, kann 

‚sowohl zum Selbststudium, das durch die zahlreichen 
guten Abbildungen und eine Formelzusammenstellung 
erleiehtert wird, wie auch zum Gebrauch neben Vor- 
lesungen empfohlen werden. 


Dresden. 
r 

Dr. Heinrich Blasius, Wärmelehre. Physi- 
kalische Grundlagen vom techni- 
schen Standpunkt. (5. Aufl.) 294 S. m. 126 
A’bb., 67 Aufg. u. 15 Tabellen. Hamburg 1949. Boysen 
& Maasch Verlag. Preis brosch. 11,50 DM. 

Die jetzt erschienene fünfte Auflage dieses Buches 
istgegenüber der vierten, die in dieser Zeitschrift Bd. 22 
(1942) S. 62 und 240 von Prof. Dr. Pa uer besprochen 
wurde, kaum verändert. Ein neues Referat erübrigt sich 
daher. 


Dresden. 


Willers. 


Willers. 


* Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen 


# 3) A 4 s 
r Verlägsanstalt (Notdruck). Preis 


Die Weiterentwieklung der Gauß-Mongeschen ER 
Differential-Geometrie im 19. und 20. Jahrhundert ist 
durch die Namen Plücker, F. Klein, Möbius, 


RENTEN PARSE ERSTEN 


% U a b 


192 Eingegangene Bi 


Beider Schriftleitung sind folgende Bücher eingegangen ( 


Instituts für Strömungsforschung und . a. d. Uniy. 

tingen), Konforme Abbildung, 359 8. 

t 276 Abb. Berlin-Göttingen-Heidelberg 1948. 
Springer-Verlag. Preis brosch. 36,— DM. 


A. Ziegler (Obering.), Einfache Formeln 
für Träger über 2—10 Felder mit be- 
liebig großen Stützweiten. MS. Bonn 
1949. Ferd. Dümmlers Verlag. Preis brosch. 7,80 DM. 


Dr.-Ing. €. Pfleiderer (Prof, a.d. Techn. Hochschule 
Braunschweig), Die Kreiselpumpen für 
Flüssigkeiten und Gase. Wasserpum- 
en, entilatoren, Turbogebläse, 
urbokompressoren. 3.Aufl. 5188, mit 


NACHRICHTEN 


Die diesjährige Tagung der Deutschen Mathematiker- 
Vereinigung findet in der mit dem 18. 9. d. J. beginnen. 
‚den Woche in Köln statt. 


In Belgrad erscheint eine neue mathematische 
Zeitschrift: Bulletin de La Sociöt& des Mathematieiens 
et Physiciens de la R.P, de Serbie, die von Herrn 


ZUSCHRIFT AN DEN HERAUSGEBER 


Zu P. Riebesell: Kritische Betrachtungen zur 
sogenannten Großzahlforschung in der Tech- 
nik und zur Anwendung mathematisch-sta- 
tistischer Methoden in der Biologie und Medi- 
zin. Z. angew. Math. Mech. 28, 226234. 


Im Heft Nr. 7/8 des 28. Bandes dieser Zeitschrift hat 
Herr Riebesell Kritik an der sog. Großzahlfor- 
schung geübt. 

Da diese in der Technik aber zu bemerkenswerten 
Erfolgen geführt hat, trifft ihre Kritik, obwohl mathe- 
matisch einwandfrei, offenbar nicht das in der Praxis 
Wesentliche. 

Zunächstist klar, daß man eine beliebige, stetige und 
eindeutige Verteilungsfunktion y = f(x) eines Merk- 
mals x durch eine Transformation x = x(£) in eine 
Normalverteilung y = e”®? des neuen Merkmals £ über- 
führen kann. Ist nämlich x = «(y) die zu f(#) inverse 
Funktion, so hat man nur x(&) =u (e7#?) anzusetzen, 
um y = e”®?zuerhalten. Die Transformation des Merk- 
mals darf also nicht damit begründet werden, daß da- 
durch eine Normalverteilung (G a u ß -Kurve) des je- 
weils vorliegenden Kollektivs erreicht werden soll, son- 
dern die Wahl des Merkmals hat so zu erfolgen, daß sie 
dem Zweck der statistischen Untersuchung entspricht. 

Untersucht man z.B. das Verhalten von Güter- 
wagen nach dem Herabrollen vom Ablaufberg eines 
Rangierbahnhofs, so wird man, wenn dabei die 
Zeitdauer von Interesse ist, welche die ablaufenden 
Wagen zum Passieren einer Weichenzone brauchen, 
die Wagengeschwindigkeit als Merkmal wählen, Inter- 
essiert jedoch die Konstruktion von Gleisbremsen, so 
ist die kinetische Energie der Wagen das zweckmäßige 
Merkmal. Ob sich bei dem einen oder dem anderen 
Merkmal evtl. eine Normalverteilung ergibt, ist bei 
dessen Wahl nicht maßgebend. 

In der Technik spielt die Normalverteilung überhaupt 
nicht die wichtige Rolle, wie in anderen Gebieten der 


"EINGEGANGENE BÜCH 


U um Dr.-Ing, e. h. Werner Schmei 


Univ. » gr 
Wi n.] 2 
mat . : ll 4 be ee er 
4. Aufl. 5828. mit 14Abb. Berlin und Heidelberg 
1947. Springer Preis brosch. 39,60 Ä 


re 


d. Techn. Hoch Berlin), Vortri 
ER Dee Matrize 
itAn ndungen 8 ndT 
Wil VIN]1558. mit 3 Abb, Berlin 1940. Akad 
Verlag. Preis brösch. 10,— DM. 5 a) 
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J. Karamata of les 
ersten Jahrganges (1 liegt vor. _ Bu 


In den letzten Jahren starben von in Deutschland 

bekannten Mathematikern: ut 
Prof. Dr. T. CGarleman, Djursholm, 
Prof. Dr. M. Bauer, Budapest. 


Statistik, wo es auf die Kennzeichnung eines Kollektivs 
durch möglichst wenige Parmeter rs Das auf- 
gezeichnete Bild der Streuungskurve o dem Tech- 
niker mehr, als ihre analytische Formel. Zeigt seine 
Kurve annähernd die Form der G a u ß-Kurve, so ist 
das dem Praktiker ein di i Zeichen etwa für 
den guten Ablauf eines ’ 
Zeigt die Verteilungskurve aber zwei (oder mehrere) 
Maxima, so hat die Aufgabe, die Kurve als Superposi- 
tion von Normalverteilungen darzustellen, keine ein- 
deutige Lösung, wie Herr R. richtig bemerkt, Auch 
sein Einwand, daß erst zu untersuchen sei, ob die 
Maxima wirklich echte Buckel sind oder innerhalb der 
Streuung eines endlichen Kollektivs i 
vom mathematischen Standpunkt aus sti 
der technischen Praxis sind derarti 
wieder nur ein diagnostischer Hinweis auf eine ge; ? 
liche Unvollkommenheit des Ferti X 
die Buckel echt sind oder nicht, kommt bei einer dann 
vorzunehmenden Wiederholung der Untersuchung 
von selbst heraus. Sind sie nun echt, so wird 
forscht, wodurch sie entstehen, und der Arbei 
solange geändert, bis sie verschwunden sind die 
angestrebte Normalverteilung möglichst erreicht ist. 
Der Standpunk des technischen Betriebsmannes 
und des Wahrscheinlichkeitstheoretikers sind in der 
Statistik also recht verschieden; aber jeder in seinem 
Bereich berechtigt. In dem Buche: ‚‚Praktische Groß- 
zahlforschung‘‘ von Dr.-Ing. Karl Daeves, Berlin, 
VDI-Verlag, 1933 findet man eine erstaunliche Fülle 
von Beispielen für die Anwendung der Großzahlfor- 
schung in der Technik, wobei der erreichte Erfolg die 
Primitivität der mathematischen Hilfsmittel recht- 
fertigt, j 
Hannover, H. v.Sanden, 


Die Zuschrift hat Herrn Prof. Dr. Riebesell 


vorgelegen. 
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